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PREFACE 


The  investigation  reported  herein  was  conducted  by  personnel  of 
“the  Soils  and  Pavements  Laboratory  (S&PL),,  U.  S,  Army  Engineer  Waterways 
Experiment  Station  (WES) ,_  as  a part  of  DA  Project  No.  UA161IQ2B52E, 
"Fragment  £ind  Projectile  Penetracion  Resistance  of  Soils,"  Task  04,  Work 
Unit  013. 

The  research  was  accomplished  by  Mr.  R.  S.  Bernard  (WES)  and 
Professor  S.  V.  Hanagud,  Georgia  Institute  of  Technology,  during  the 
period  July  1973  to  June  197*+  under  the  supervision  of  Dr.  J.  G. 

Jackson,  Jr.,  Chief  of  the  Soil  Dynamics  Division,  S&PL.  Messrs.  J.  P. 
Sale  and  R.  G.  Ahlvin  were  Chief  and  Assistant  Chief,  S&PL,  respec- 
tively. SPU  D.  C.  Creighton  executed  the  necessary  cpmputer  calcula- 
tions, and  Dr.  B.  Rohani  provided  continuing  support  through  his  famil- 
iarity with,  previous  work  in  this  area.  This  report  was  prepared  by 
Mr.  Bernard  and  Professor  Hanagud. 

BG  E.  D.  Peixotto,  CE,  and  COL  G.  H.  Hilt,  CE,  were  Directors  of 
WES  during  this  investigation  and  the  preparation  of  the  report. 

Mr.  F.  R.  Brown  was  Technical  Director. 
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COHVERSIONi  FACTORS,,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  Sv  customary  units,  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Si)'  units  as  follows: 


Multinl 


To. Obtain 


inches 

2.54 

centimetres 

feet 

0.3048 

metres 

pounds  (mass) 

0.4535924 

kilograms 

slug-square  feet 

1.353 

kilogram-square  metres 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  metre 

slugs  per  cubic  foot 

515.3788 

kilograms  per  cubic  metre 

pounds  (force)  per  square  inch 

6894. 75T 

pascals 

kips  (force),  per  square  inch 

6894.757 

kilopascals 

feet  per  second 

0.3048 

metres  per  second 

degrees  (angle). 

0.0174539 

radians 

Note:  1 kbar  = 1,000  bars  = 14,500  pounds  per  square  inch 

= 100  megapascals. 

1 slug  = 32.2  pounds  (mass). 
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-DEVELOPMENT  OF  A PROJECTILE  PENETRATION  THEORY 
PENETRATION  THEORY  FOR  SHALLOW  TO  MODERATE  DEPTHS 

PART  I'!  INTRODUCTION 
Background 

1.  The  phenomenon  of  projectile  penetration  has  been  the  object 
of  numerous  analytical  and  experimental  investigations  over  the  past 
two  centuries. Until  recent  years,  successful  analyses  were  of  an  em- 
pirical or  semiempirical  nature,  utilizing  various  target  resistance 
functions  and  pseudoconstants  which  were  inevitably  drawn  directly  from 
penetration  experiments.  As  a result,  the  subject  of  penetration  re- 
mained isolated  from  the  rest  o^  mechanics  because  the  relation  between 
target  constitutive  properties  and  penetration  resistance  was  unknown. 
Nevertheless,  a number  of;  useful  empirical  equations  were  developed 

which  afforded  reasonable  ..estimates  of  penetration  depth  for  a variety 
1-4 

of  situations.  Unfortunately,  the  reliability  of  such  equations  is 
limited  to  the  range  of  test  conditions  for  which  they  were  deduced,  and 
the  accuracy  of  prediction  may  dwindle  significantly  with  extrapolation 
to  situations  which  have  not  been  characterized  experimentally. 

2.  The  last  decade  has  given  rise  to  a more  basic  approach  to 

penetration  through  the  use  of  digital  computers.  In  various  efforts, 

the  details  of  target  and  projectile  motion  have  been  analyzed  by  means 

of  two-dimensional  finite-difference  techniques  employing  realistic  con- 

5-7 

stitutive  properties  for  both  target  and  projectile.  While  a finite- 
difference  solution  may  achieve  good  agreement  with  experimental  observa- 
tion, as  shown  in  Reference  7S  its  execution  is  cumbersome  and  costly, 
thus  diminishing  its  utility  as  a practical  engineering  tool.  However, 
this  approach  represents  a powerful  device  for  examining  the  details  of 
target  and  projectile  behavior,  indicating  the  relative  consequences  of 
different  simplifying  assumptions,  and  providing  a means  of  detei’mining 
what  mechanisms  are  most  important  in  the  projectile-target  interaction. 
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3.  At  the  present 'time,  there  Is  a growing  heed  for  a simple 
theory  which  affords  reliable  analyses  of  the  impact  and  penetration  of 
rigid  projectiles  such  as  missiles,  bombs,  bullets,  shell,  fragments, 
artillery  and  mortar  rounds  * and  air-delivered  mines  and  sensors  into 
targets  such  as  soil,  rock,  concrete,  and  metal.  The  results  Of  these 
analyses  could  be  used  in  design  and  optimization,  both  from  the  offen- 
sive and  defensive  points  of  view.  In  July  1973,  a small  research  ef- 
fort was  initiated  under  the  sponsorship  of  the  Office,  Chief  of  Engi- 
neers,, to  develop  a,  tractable  projectile  penetration  theory  fOr  natural 
and  artificial  earth  targets.  The  general  requirements  were  that  the 
theory  should  account  for  the  penetration  of  homogeneous  and  layered 
targets  by  rigid  projectiles  of  various  nose  shapes  and  should  deal  at 
least  approximately  with  oblique  impact.  Further,  the  theory  should  be 
derived  from  a consideration  of  target  constitutive  properties  without 
recourse  to  penetration  tests  per  se.  The  objective  was  to  replace- 
neither  empirical  nor  finite-difference -analyses,  but  rather  to  provide 
a practical  tool-  for  making  reasonable  predictions  before  the  execution 
Of  a penetration  test  or  finite-difference  calculation.  The  develop- 
ments presented  in  this  report  are  the  result  of  efforts  conducted  under 
this  program  during  the  period  July  1973  to  June  197**.  Report  2 will 
document  the  extension  of  the  present  theory  for  application  to  very 
deep  penetration  in  homogeneous  and  layered  targets. 

Purpose 

4.  The  purpose  of  this  investigation  is  to  build  upon  existing 

8-12 

penetration  theory  to  develop  an  approximate  but  more  general  theory 
meeting-  the  requirements  stated  above.  The  starting  point  for  the  pres- 
ent work  is  the  so-called  "cavity  expansion  theory,"  originally  applied 

9 10 

to  projectile  penetration  by  Goodier. 


5.  Part  II  of  this  report  outlines  modifications  of  previous 
8-12 

theoretical  work  ‘ in  which  the  effects  of  different  nose  shapes, 
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layered  targets,  and  oblique  impact,  are  incorporated  into  a penetration 
theory  for  rigid  projectiles.  In  Part  III,  the  modified,  theory  is  non- 
dimensionalized,  governing  parameters  are  determined,  and  results  of 
parameter  studies  are  presented.  Part  IV  contains  comparisons  of  theo- 
retical results'  with  experimental  observations  and  two-dimensional 
finite-difference  solutions.  Part  V presents  the  conclusions  and  recom- 
mendations, The  theory  of  dynamic  expansion  of  a spherical  cavity  in  a 
concentrically  layered  medium  is  developed  in  Appendix  A.  A discussion 
of  material  compressibility  as  it  relates  to  the  cavity  expansion  theory 
and  the  penetration  theory  is  contained  in  Appendix  B.  The  effect  of 
target  acceleration  on  final  penetration  depth  predictions  is  investi- 
gated in  Appendix  C. 
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PART  II:  PENETRATION  THEORY 

Historical  Background 

6.  The  foundation  for  the  present  penetration  theory  was  laid  by 

Q 

Bishop,  Hill,  and  Mott,  who  applied  the  theory  of  spherical  cavity  ex- 
pansion in  a strain-hardening  medium  to  penetration  by  a rigid  static 
punch.  Goodier^  later  extended  this  work  and  applied  it  to  the  pene- 
tration of  a strain-hardening  target  by  a high-velocity  rigid  spherical 

12 

projectile.  Hanagud  and  Ross  modified  the  cavity  expansion  theory  to 
account  for  a compressible  medium  by  means  of  a locking  approximation 
and  employed  the  results  in  Goodier's  penetration  theory. 

Cavity  Expansion  Theory 

7.  The  theory  of  dynamic  expansion  of  a spherical  cavity  sur- 
rounded by  concentric  layers  of  compressible  locking,  strain-hardening 
material  is  developed  in  Appendix  A.  This  cavity  expansion  theory  forms 
the  basis  for  the  penetration  theory  discussed  in  the  present  report. 

The  behavior  of  the  material  in  shear  (Figure  A2)  is  characterized  by 
the  compressive  yield  strength  Y , the  modulus  of  elasticity  E , and 
the  strain-hardening  modulus  E^  .*  The  effect  of  material  compressi- 
bility is  approximated  by  means  of  an  instantaneous  locking  assumption 
whereby  material  in  the  plastic  state  is  presumed  to  be  compressed  to  a 
uniform  (i.e.  locked)  density  . As  discussed  in  Appendix  B,  the 
instantaneous  locked  density  p^  is  a function  of  the  instantaneous 
average  dynamic  pressure  in  the  plastic  zone. 

8.  According  to  Equation  A76,  the  compressive  normal  stress  p 
at  the  cavity  surface  is 

p = ps  + P*(Biaa  + S ps  + PI  ^ 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de- 
fined in  the  Notation  (Appendix  D). 
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The  two  terms,  p and  pT  , which  appear  in  Equation  1 are  the  sepa- 

Si' 

rate  contributions  due  to  material  deformation  (shear)  and  inertia,,  re- 


spectively. Henceforth,  the  shear  contribution  p 

s 

"shear  resistance";  and  the  inertial  contribution  p 


will  be  called  the 

i = "A" ♦ V2) 


will  be  called  the  "dynamic  pressure"  where  a , a , and  S are  the 


radial  position,  velocity,  and  acceleration  of  the  cavity  wall,  and 
and  B2  are  dimensionless  inertial  coefficients.  For  a homogeneous  ma- 
terial (i.e.  no  layering),  p , B..  , and  B0  are  identified- by  Equa- 

tions  A87,  A88,  and  A89,  respectively.  Moreover,  these  quantities  are 


functions  of  Y , E , , and  the  instantaneous  locking  strain 

= In  p/p^  » where  p is  the  initial  density. 


Penetration  Theory  for  Spherical  Projectiles 


9.  While  spherical  cavity  expansion  and  projectile  penetration 
can  hardly  be  regarded  as  the  same  process,  the  penetration  of  a solid 
target  by  a rigid  spherical  projectile  (Figure  l)  suggests  at  least  a 
crude  analogy  with  spherical  cavity  expansion.  Gooaier  is  the  first 
investigator  to  apply  such  an  analogy  to  projectile  penetration.  Con- 
sidering only  rigid  spherical  projectiles,  he  makes  the  following 
assumptions : 
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a.  The  embedded  portion  of  the  projectile  frontal  surface  as 
seen  in  the  above  figure  is  in  complete  contact  with  the 
target . 

b.  During  the  embedding  process  (penetration  depth  q <_ 

projectile  radius  aQ),  the  projectile-target  'nteraction 

is  equivalent  to  a dynamic  Brinell  hardness  test*  and  the 

9-12 

s al  resisting  force  is  given  by  Meyer's  law. 
c_.  After  the  projectile  is  fully  embedded  (q  > a^),  only 
normal  stresses  exist  on  the  projectile  frontal  surface 
(i.e.,  the  effect  of  friction  is  negligible). 

10.  At  this  point,  Goodier  sets  about  to  modify  Equation  1 in 
order  to  approximate  the  normal  stress  distribution  on  the  frontal  sur- 
face of  a fully  embedded  spherical  projectile.  To  this  end  he  makes  two 
assumptions:  (a)  the  normal  stress  due  to  shear  (i.e.  shear  resistance) 

is  constant  over  the  frontal  surface  and  is  identical  with  pg  in  Equa- 
i tion  1;  and  (b)  p^  is  maximum  at  the  stagnation  point  0=0  and  zero 

■ at  the  shoulder  0 = ir/2  , where  0 is  measured  as  shown  in  Figure  1. 

; Identifying  aQ  with  the  cavity  radius  a , the  projectile  velocity 

q with  the  cavity  radial  velocity  a , and  the  projectile  acceleration 
| § with  the  cavity  radial  acceleration  £ , Goodier  then  proposes 

PI  = Pt(BlV*  + B2^2)  COS  9 

as  a reasonable  variation  of  dynamic  pressure  over  the  frontal  surface 
0 £ 0 < ir/2  . The  compressive  normal  stress  p on  the  projectile 
frontal  surface  is  then 

| 

•;  p * ps  + ^(Vo5  + M2)  005  0 (3) 

j>.  In  Goodier1  s investigation  the  target  is  assumed  incompressible  such 

, that^’10  p0  = p , B = 1 , and  B«  = 3/2  . The  rationale  employed  by 

‘ ' ’ * 11  12-1  d 
< Hanagud  and  Ross  ’ is  identical  with  that  of  Goodier  except  that  the 

V , target  is  considered  compressible  locking  such  that  p^  > p , B^  < 1 , 

j B0  < 3/2  , and  the  value  of  p is  less  than  the  incompressible  value. 
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11.  Integration  of  Equation  3 over  the  frontal  surface  leads  to 
the  following  equation  of  motion  for  a projectile  of  mass  M : 


Mq  = -*«^PS  + § Pt(V05  + B2«2)] 


Equation  1 is  rearranged  such  that 


(M  + I "4’ A)*  ■ -*ao(l>s  * I C4B2* 


Equation  A88  indicates  that  B <_  1 ; consequently,  the  effect  of  target 

13 

acceleration  is  small  in  Equations  4 and  5 when  the  added  mass'  for  the 
sphere  is  small  compared  to  the  projectile  mass,  i„e. , when 


| « M 


With  the  ratio  R defined  as 
a 


Vt 


the  effect  of  target  acceleration  on  projectile  loading  is  then  small  v 
when  R « 1 . Moreover,  the  quantity  (2/3)ita  p.B1  probably  repre- 
sents  a high  upper  limit  on  the  added  mass  term,  and  penetration  pre- 
dictions are  improved  in  some  cases  by  discarding  the  added  mass  term 
(i.e.  the  term)  altogether,  as  shown  in  Appendix  C.  At  any  rate, 
Goodier’s  dynamic  pressure  distribution  is  adequately  approximated  by 


PI  ~ p£B2q  C0S  9 » Ra  <<:  1 


12.  Although  the  rationale  employed  by  Goodier  leads  to  fairly 
accurate  predictions  of  final  penetration  depth  (Appendix  C),  the 
cos  0 distribution  of  p^  is  somewhat  arbitrary.  Other  distribu- 
tions of  dynamic  pressure  produce  projectile  equations  of  motion  iden- 
tical to  Equations  4 and  5-  (For  example,  this  occurs  when  cos  .6  is 

4 \ 

replaced  by  1 - sin  0 in  Equations  2 and  3.)  Furthermore,  the 
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projectile  equation  of  motion  must  be  integrated  with  respect  to  time  to 
obtain  the  final  depth  of  penetration,  so  that  a "correct”  prediction  of 
final  depth  alone  indicates  only  that  the  average  resisting  force  is 
correctly  predicted.  Thus,  there  is  some  flexibility  in  Goodier's 
rationale,  insofar  as  the  distribution  of  dynamic  pressure  is  concerned. 
Moreover,  it  is  this  flexibility  which  allows  extension  of  the  cavity 
expansion  analogy  to  projectiles  with  nose  shapes  other  than 
hemispherical. 
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Penetration  Theory  for  Axisymmetric  Projectiles 

13.  Before  generalizing  the  penetration  theory  to  nose  shapes 
other  than  hemispherical,  it  is  necessary  to  reevaluate  Goodier's  three 
basic  assumptions.  First,  the  assumption  of  complete  contact  between 
the  target  and  the  embedded  portion  of  the  projectile  nose  (i.e.  frontal 
surface)  is  retained.  Although  this  assumption  is  physically  question- 
able, it  will  provide  an  upper  limit  on  the  resisting  force,  at  least 
within  the  framework  of  the  theory.  Second,  the  dynamic  Brinell  hard- 
ness analogy  is  discarded.  Instead,  an  analogy  with  the  cavity  expan- 
sion theory  will  be  used  throughout  the  penetration  process  to  approxi- 

f 

mate  the  normal  stress  distribution  on  the  embedded  portion  of  the 
projectile  nose.  (This  will  produce  a higher  target  resistance  during 
the  embedding  process,  which  will  improve  penetration  depth  predictions, 
especially  for  shallow  penetration.)  Third,  friction  effects  are  dis- 
regarded at  all  times  on  the  projectile  nose  and  aft  body  alike.  This 

is  a questionable  assumption,  although  there  is  some  evidence  suggesting 

7 lU 

that  friction  on  the  aft  body  is  relatively  insignificant. 

14.  Subject  to  the  above  assumptions.  Equation  1 is  now  modified 
in  order  to  approximate  the  normal  stress  distribution  on  the  nose  of  a 
projectile  with  an  arbitrary  axisymmetric  nose  shape.  As  in  Goodier's 
rationale,  the  normal  stress  due  to  shear  is  assumed  uniform  over  the 
nose  and  equal  to  p . This  leaves  only  the  distribution  pT  to  be 
determined  as  a function  of  nose  shape.  The  problem,  of  course,  is  to 
arrive  at  a single  rationale  which  generalizes  the  distribution  of  Pj 
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to  any  axisymmetric  nose  shape  in  a straightforward  manner-  Further- 
more, the  resulting  effect  of  nose  shape  on  final  penetration  depth 
should  he  consistent  with  .experimental  observation. 

15..  One  approach  is  simply  to  guess  the  form  of  the  pT  dlstri- 

o 1q  A 

hution  in  a manner  similar  to  Goodier.  * Consider  a projectile  with  a 
fully  embedded  conical  nose,  as  shown  in  Figure  2.  „ As  in  the  case  of  a 


Figure  2.  Projectile  with  a conical  nose 

hemispherical  nose,  one  might  assume  that  p^.  = Pfc(Biaoq.  + B^q2^  at  the 
stagnation  point  (i.e.,  at  the  tip  of  the  nose,  z * 0).  At  the  base  of 
the  nose,  z * L , one  might  assume  that  the  dynamic  pressure  reduces  to 


*2\  xi 

B^aQq  + Bgq  J Bin  0 . This  value  corresponds  to  the  value  of 
Pj  cn  a hemispherical  nose  at  polar  angle  0 = ir/2  - 0 when  p^.  is 
distributed  according  to  cos  0 (as  opposed  to  cos  0).  The  variation 
of  Pj  from  z * 0 to  z = L could  then  be  expressed  according  to 

pT  = p^B^q  + B2q2^£sinn  0 + (l  - sin11  0)f(zjj  (9) 


where  f(z)  is  a dimensionless  function  which  decreases  smoothly  from 
unity  at  z = 0 to  zero  at  z = L . 

16.  Since  the  problem  at  hand  is  to  make  an  ad  hoc  modification 
of  the  dynamic  pressure  for  different  nose  shapes,  there  is  basically 
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nothing  wrong  with  the  approach  outlined  above.  On  the  other  hand,  the 
arbitrariness  in  the  development  of  Equation  9 prompts  one  to  seek  a 
more  fundamental  scheme  based  (albeit  crudely),  on  the  underlying  physics. 
The  fact  that  the  dynamic  pressure  which  appears  in  Equation  1 depends 
on  both  particle  acceleration  and  velocity  suggests  that  the  dis- 
tribution for  an  axisymmetric  projectile  can  be  obtained  by  approxi- 
mating the  variation  of  particle  acceleration  and  velocity  along  the  sur- 
face of  the  nose.  Moreover,  since  most  cases  of  practical  interest 
satisfy  the  condition  R « 1 , the  particle  acceleration  has  only  a 
small  effect  on  projectile  loading,  and  it  is  then  sufficient  to  ap- 
proximate the  variation  of  particle  velocity  alone. 

17 . In  order  to  construct  a reasonable  distribution  of  particle 
velocity  v on  the  nose  surfac^,  one  must  first  consider  certain 
physical  constraints  on  particle  motion: 

a.  A particle  at  the  very  tip  of  the  nose  (i.e.  the  stagna- 
tion point)  must  move  with  the  projectile  velocity  q • 

b.  The  component  of  particle  velocity  v normal  to  the 
surface  of  the  nose  must  be  equal  to  the  normal  component 
of  projectile  velocity  (i.e.,  material  cannot  cross  the 
target/projectile  interface). 

c.  Over  the  nose  surface  v must  be  continuous. 

P 

d.  Finally,  it  is  assumed  that  the  tangential  component  of 
particle  velocity  v^  is  zero  at  the  base  of  the  nose. 

Relative  to  these  constraints,  only  the  variation  of  v^  along  the 
length  of  the  nose  remains  to  be  specified.  Although  it  is  quite  pos- 
sible that  Vj.  is  a complicated  function  of  target  properties,  projec- 
tile velocity,  nose  shape,  and  relative  position  on  the  nose,  the  pres- 
ent investigation  considers  only  the  effects  of  nose  shape  and  relative 
position. 

18.  Consider  first  the  fully  embedded  conical  nose  shown  in 
Figure  2.  The  normal  component  of  velocity  v^  is  given  by 

v = a sin  0 
n ^ 

' Ik 


(10) 


where  f is  the  half-angle  for  the  circumscribed  cone  vMch  is  tangent 
to  the  ogive  at  any  longitudinal  position  z . To  arrive  at  a func- 
tional form  for  v^  which  is  consistent  with  Equation  12,  the  effective 
cone  length  L (f)  and  effective  longitudinal  position  z ($)  (Fig- 
ure  3)  are  defined  by 


- cos  0 cos 


sin 


where 


z(i>)  H a(^)  cot  ^ 


a(\J»)  - R(cos  - cos  0) 


a = R(l  - cos 
o 


The  effective  dimensionless  cone  position  x (^)  is  then  defined  as 


L (ij/)  1 - cos  0 cos 


In  order  to  be  consistent  with  Equation  12,  is  taken  to  be 


= qjl  - xe(*)]1/2  cos  i|» 


Thus,  the  expression  -for  is 


(vn  + vt)  * " xe^  cos2 


20.  By  replacing  a by  aQ  , a by  v^  , and  a by  q in 
Equation  1,  the  general  expression  for  p,.  on  an  axi symmetric  nose  sur 
face  is 


PI  “ W + plVp  “ PlB2Vp  ’ Ra  <<  i 


where  v^  is  given  by  Equations  13  and  21  for.  cones  and  ogives,  respec- 
tively. Thus, 

2 *2 

Pj  ® p^BgCl  - x cos  0)q  for  cones  (23) 


PI  * cos2  *]42  for  ogives  (2U) 

vhen  R « 1 . For  a hemispherical  nose,  which  is  the  limiting  case 
for  an  ogive.  Equation  24  reduces  to 

PI  ~ PtB2^2^  “ CO®  ^ “ P£B2q2(l  - sin*  6)  (25) 

where  0 - ir/2  - iji  as  indicated  in  Figure  1.  A comparison  of  Equa- 
tions 8 and  25  is  presented  in  Figure  4 to  show  that  for  a hemispherical 
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v — PRESENT  WORK,  EQUATION  25 
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Figure  4.  Graphic  comparison  of  Equations  8 and  25  for  distribution 
of  dynamic  pressure  oh  a hemispherical  nose 


nose  the  p^,  distribution  obtained  by  assuming  a functional  form  for 
Vj.  is  not  radically  different  from  Goodier's  assumed  Pj  distribution. 

21.  With  the  normal  stress  distribution  having  been  approximated 
for  conical  and  ogival  noses  (the  rationale  outlined  in  paragraphs  17-20 
can  be  readily  applied  to  other  axi symmetric  nose  shapes) , the  next  step 
is  to,  obtain  the  resulting  projectile  equations,  of  motion.  This  in- 
volves integrating  the  normal  stress  over  the  surface  of  the  nose  to  ob- 
tain the  total  axial  resisting  force  F acting  on  the  projectile. 

z 

Thus,  for  projectiles  with  conical  noses.  Equations  11,  13,  and.  22  are 
substituted  into  Equation  1 and  integrated  over  the  surface  of  the  nose 
to  obtain 

* z=L 

Mq  - -7  z~  - J [ps  + PI(z»0)]J;  sin  0 dA  (z,0);  (26) 

z=0 

where 

dA  (z,0)  = 2rr  ■§•  z dz  (27) 

cos  0 

Likewise,  for  projectiles  with  ogival  noses.  Equations  19,  21,  and  22  are 
substituted  into  Equation  1 and  integrated  over  the  surface  of  the  nose 
to  obtain  , 


where 

dA  (f)  = 2rr  &£JL  z (,),)  dz  (^)  (29) 

cos  ij/  e e 


(i(>)J  sin  tj>  dA  (ip)  (28) 


The  resulting  projectile  equations  of  motion  have  the  form 

(M  + iraVlOq  ='-ira?fp„  + paBJ  q2)  (30) 


The  quantity  (no3e  shape  factor)  results  from  the  variation  of  p^_ 
over  the  nose  and  thus  characterizes  the  relative  effect  of  nose  shape 
on'  the  average  inertial  resistance'  of  the  target  to  penetration.  For 
conical  noses  of  length  L and  base  diameter  D = 2aQ  shown  in  Fig- 
ure 2,.  f is  expressed  as 


f 


n 


Is! 


+ ly 


(31) 


where 


C = | ‘ (32) 

For  ogival  noses  of  length  L and  base  diameter  D illustrated  in  Fig- 
ure 3,  f is  written  as 

f = 1 r g-  [b2  In  (2  - e)  - (3B2  + 2B)(B  - e) 

e (1  - e)°  L 

+ | (3B2  + 6B  + 1)(B2  - e2)  - | (B5  - e5) 

- | (B2  + 6B  + 3)(B3  - e3)  + (2B  + 3)(B^  - e*4)]  (33) 

where 

B = 2e  - e2  (34) 

e 5 (35) 

4c2  + l 

and  ? = L/D  . For  ogives,  the  relation  between  the  caliber  radius  head 
(CRH)  and  x,  is 


2 1 
CRH  * 5 + f 


(36) 


An  ogive  with  C = 1/2  is  a hemisphere,  for  which  case  Equation  33  re- 
duces to  f = 2/3  and  the  right  hand  side  of  Equation  30  is  identical 


with  Equation  5.  Thus,  for  spherical  projectiles  and  projectiles  with 

hemispherical  nose  shapes,  the  distribution  of  particle  velocity  given 

by  Equation  21  leads  to  essentially  the  same  equation  of  motion  as  that 

q in  ji  ip 

obtained  by  Goodier  * and  by  Hanagud  and  Ross  * when  R « 1. 

& 

22.  Whether  one  chooses  to  develop  projectile  equations  of  motion 

based  on  assumed  stress  distributions  or  assumed  particle  velocity  dis- 
tributions is  of  little  consequence  in  the  final  analysis.  The  impor- 
tant question  is  whether  or  not  the  equation  of  motion  (30)  leads  to 
reliable  predictions  of  final  penetration  depth  for  projectiles  with 
different  nose  shapes.  Moreover,  the  relative  effect  of  nose  shape  on 
final  penetration  depth  should  lie  within  experimentally  observed  bounds. 
(For  earth  penetration,  variation  of  nose  shape  from  a blunt  cylinder  to 
a sharp  cone  seldom  increases  final  penetration  depth  by  more  than  a 
factor  of  three.  ) Thus,  the  choice  of  a velocity  or  stress  distribu- 
tion leading  to  f is  neither  arbitrary  nor  unique,  and  it  is  no  acci- 

dent that  variation  of  nose  shape  from  a blunt  cylinder  to  a sharp  cone 
reduces  target  inertial  resistance  in  Equation  30  by  no  more  than  a fac- 
tor of  three.  The  tangential  velocity  distribution  given  by  Equation  12 
is  chosen  so  as  to  produce  precisely  this  effect. 

23.  While  f is  presumed  to  characterize  the  relative  effect  of 

n 

nose  shape  on  inertial  resistance  to  penetration  for  a fully  embedded 
nose,  the  embedding  process  itself  must  also  be  considered  in  order  to 
determine  the  overall  effect  of  nose  shape  on  penetration  depth,  espe- 
cially for  shallow  penetration.  Assuming  that  f approximately  char- 
acterizes inertial  resistance  during  the  embedding  of  the  nose,  the 
equation  of  motion  during  the. embedding  process  is 

[m  + TO^(q)pABjq  = -ira^(q)|ps  + p£B2fn^2]  » <1  £ L (37) 

where  a (q)  is  the  cylindrical  radius  of  the  nose  at  the  surface  of 
® 2 

the  target,  and  rar(q)  is  the  frontal  area  of  the  projectile  which  is 
in  contact  with  the  target  for  a given  penetration  depth  q <_  L , as 
shown  in  Figure  5.  When  q > L , Equation  37  is  replaced  by  Equation  30. 
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Figure  5.  Projectile  with  partially  embedded  nose 


A study  of  the  relative  effect  of  nose  shape  on  final  penetration  depth 
is  given  in  Part  III, 

Penetration  Theory  for  Layered  Targets 

24.  The  penetration  theory  for  layered  targets  is  based  upon  the 
cavity  expansion  theory  for  a concentrically  layered  medium  presented  in 
Appendix  A.  The  physical  situation  to  be  considered  is  illustrated  in 
Figure  6.  As  is  the  case' in  Appendix  A,  the  investigation  is  restricted 
to  the  two-layer  situation.  Physically,  this  means  that  a target  may 
have  any  number  of  layers  so  long  as  the  projectile  equation  of  motion 
is  significantly  influenced  by  no  more  than  two  layers  at  a time. 

25.  The  projectile  equation  of  motion  is  assumed  to  have  the 
same  form  as  that  for  a homogeneous  target.  Equation  30.  However,  ap- 
proximations for  pg  , B1  , and  Bg  must  be  obtained  through  an  appro- 
priate analogy  with  the  cavity  expansion  theory  for  a concentrically 
layered  medium.  The  point  which  becomes  the  point  of  maximum  deforma- 
tion on  the  layer  interface  is  denoted  by  H(t)  and  coincides  with  the 
projectile  axis  of  symmetry  (Figure  6).  The  initial  distance  from  the 


Figure  6.  Penetrabion  of  a layered  target 

target  surface  to  the  layer  interface  is  H.  . A first  approximation 

o 

(probably  representing  an  upper  limit)  to  the  effective  deformation  of 
the  layer  interface  is  then  obtained  by  replacing  the  penetrating  pro- 
jectile with  a spherical  cavity  of  radius  aQ  and  radial  velocity 
4 located  at  the  base  of  the  projectile  nose  (i.e.,  at  depth  q - L) 
and  assuming  that  the  motion  of  the  point  H(t)  is  the  same  as  that 
of  a point  lying  on  a spherical  surface  of  instantaneous  effective 
radius  H - q + L and  effective  radial  velocity  H , concentric  to 
the  effective  cavity  (a  = aQ).  The  motion  of  the  plastic  front  S(t) 

can  be  determined  in  the  same  way,  where  S - q + L is  the  instanta- 

' « 

neous  effective  radius  of  the  front  and  S is  its  effective  radial 
velocity.  First  approximations  to  p , B.  , and  B_  are  then 

3 X ct 
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obtained  by  substituting  aQ  for  a(t)  , q for  a , - q + L for 

h , H - q + L for  h{t)  , H for  h , S - q + L for  b(t)  , and 
S for  b in  Equations  A79-A8U . 

26.  In  direct  analogy  with  Appendix  A,  H s Hq  for  S <_  Hq  , 
i.e.,  the  deformation  of  the  interface  is  assumed  negligible  prior  to 
the  arrival  of  the  plastic  front.  From  Equation  A57,  S - q + L is  re- 
lated to  ao  by 


S - q + L = 


,1/3 


S < H 


(38) 


where  6^  is  given  by  Equation  A53.  The  condition  S £ Hq  is  now  re- 
placed by  q £ q*  where 


> t 


9*  = H0  + L - 


*1/3 


(39) 


which  represents  projectile  penetration  depth  when  the  plastic  front 
reaches  the  layer  interface.  Equation  39  is  structured  so  that  S = 
when  q = q#  . When  q <_  q*  , the  projectile  equation  of  motion  is 


•5  (n  ) 

M + 


q = —no. 


(1)  (1)  -2 

5s  + \B2  fnq 


q < q* 


(**o) 


where,  with  the  exception  of  and  , the  subscripts  1 and 

2 denote  first  and  second  layer  properties,  and 


j(l)  = V ’’nd  + — E 
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(41) 
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B. 


(1) 


- 1 - ^/3 + f - D 


,1/3 
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'I 


H - q + L 
o 


and 


B, 
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(43) 


where  6 , 8 , and  a are  defined  by  Equations  A90,  A91,  and  A93,  re- 

spectively. Equations  4l-43  are  analogous  to  Equations  A79-A81. 

27.  After  the  plastic  front  has  reached  the  interface,  the  pro- 
jectile equation  of  motion  is 


L*  + ”aopl1Bi*)]’  = + P11B2*)V2]  • 


q > q» 


m. 


(*)  (*) 

where,  with  the  exceptions  of  B^  and  B^  , the  subscripts  1 and 

2 denote  first  and  second  layer  material  properties,  £pd 
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(45) 


24 


p„  o,(l  - 6n)‘ 


,^jt 

\S  - q + L/  / 


2Pjli  (1  - o2] 


Equations  45-47  are  analogous  to  Equations  A82-A84.  Examining  Equa 
tion  A64,  and  substituting  H - q + L for  h , H for  h , S - q 
for  b , and  S for  b , then 

(H  - q + L)2H  = 62(S  - q + L)2S  * a2q  , q > q# 


vhere  is  given  by  Equation  A6l.  Equation  48  leads  to 


H = a 


I ^ 2 + H0  , q> 

J (H  - q’  + L)2 


(S  - q’  + L) 


Equations  1+9  and  50  give  only  crude  (and  probably  upper  limit)  approxi- 
mations to  the  effective  positions  of  the  layer  interface  H and  the 
plastic  front  S as  the  projectile  approaches  the  interface.  Moreover, 

Equation  50  is  reasonable  only  when  S_  < 6.  , i.e.,  when  the  distance 

^ ^ 1/3 

S - q + L increases  to  the  final  value  a Q/6g  as  the  projectile  en- 
ters the  second  layer.  On  the  other  hand,  when  6 > 6 , the  distance 

1/3  2 1 

S - q + L decreases  to  the  final  value  a /S'  , in  which  case  S 

0 2 1/3 

is  simply  located  at  H until  S - q + L = &Q/&2  * 

28.  The  approximations  engendered  by  placing,  the  effective  cavity 
(a  - aQ)  at  the  base  of  the  nose  require  some  modification  when  q > H , 
i.e. , when  the  nose  extends  beyond  the  effective  interface  depth  H , as 
illustrated  in  Figure  7. 
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Figure  7*  Situation  when  q > H 


In  such  a situation,  the  frontal  area  in  contact  with  layer  2 is  m 


and  the  frontal  area  in  contact  with  layer  1 is  ir^a2  - a|^  . 
projectile  equation  of  motion  then  becomes  approximately 


H 


The 
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for  q > H , where 
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The  quantities  p , B , and  B'  are  shear  resistance  and 

S X u 

inertial  coefficients  for  a semi-infinite  second  layer,  given  by 
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(2)  2 1* 
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(60) 


29.  When  the  projectile  nose  becomes  essentially  embedded  in  the 
second  layer,  i.e.,  when  H - q + L - aQ  , the  equation  of  motion  is 
approximately 


3 (2)1  ••  2 r 

M + ira^p.  B.v  ; q S£  -ira  ] 
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(61) 


which  is  the  same  as  for  a homogeneous  semi-infinite  target  having  the 
same  properties  as  the  second  layer. 

30.  When  q < L , i.e.,  when  the  nose  is  not  fully  embedded  in 
the  first  layer  of  the  target,  Equations  38-6l  still  apply  except  that 
aQ  is  replaced  by  a^  and  L is  replaced  by  q.  (Figure  5). 

Oblique  Impact  Theory 

31.  The  investigation  thus  far  has  disregarded  projectile  rota- 
tion which  may  occur  due  to  oblique  incidence  with  the  target.  Conse- 
quently, the  next  step  involves  a projectile  with  a conical  nose  of  half- 
angle 0 which  strikes  the  target  with  an  initial  angle  of  obliquity 

, as  illustrated  in  Figure  8,  The  X,Y,Z  right-hand,  orthogonal 
coordinate  system  is  fixed  at  the  point  of  impact  on  the  surface  of  the 
target.  Figure  9 shows  the  x,y,z  right-hand,  orthogonal  coordinate 
system  fixed  on  the  projectile  nose  tip.  The  z axis  coincides  with  the 
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Figure  8.  Oblique  impact  of  a pro-  Figure  9*  Coordinate' systems 
Jectile  with  a conical  nose 


projectile  axis  of  symmetry,  and  at  impact  the  projectile  y axis  coin- 
cides with  the  fixed  Y axis.  With  to  as  the  instantaneous  angle  of 
obliquity,  the  projectile  may  undergo  rotation  during,  the  embedding 
process  due  to  asymmetric  loading  over  the  shaded  portion  of  the  nose 
(Figure  10).  The  rotation  of  the  projectile  about  its  center  of  gravity 
(CG)  is  now  investigated,  subject  to  the  following  restrictions: 

a.  The  investigation  is  limited  to  the  embedding  of  the  nose, 
and  projectile-target  interactions  aft  of  the  nose  are 
not  considered. 

b.  The  effect  of  angle  of  attack  upon  the  stress  distribu- 
tion on  the  projectile  is  disregarded. 

c.  The  location  z_  of  the  CG  is  aft  of  the  nose,  i.e. 

VJ 

zQ  is  greater  than  L , as  shown  in  Figure  10. 

d.  The  embedded  surface  of  the  nose  is  assumed  to  be  every- 
where in  contact  with  the  target. 

e_.  Finally,  the  compressive  normal  stress  pn  on  the  em- 
bedded surface  of  the  nose  is  assumed  to  be 

p as  p + p B_f  q2  (62) 

fs  H 2 r 

where  q is  the  translational  velocity  of  the  CG  . 
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Figure  10.  Projectile  with  conical  nose  partially  embedded  at 

angle  of  obliquity  u 

Equation  62  corresponds  approximately  to  the  average  compressive  normal 
* -v  stress  on  a partially  embedded  nose  when  R « 1 (Equation  37). 

32.  For  computation  of  forces  and  moments  acting  on  the  projec- 
tile, the  cone  is  divided  into  circular  bands  of  radius  r(z)  and  in- 
finitesimal height  dz  . However,  in  the  shaded  region  (Figure  10 ), 

only  a portion  of  each  band  is  subjected  to  p , and  for  each  band 

n 

r(z)  there  is  a limiting  azimuthal  angle  (z)  bounding  the  arc 

max 

over  which  pn  acts.  Looking  along  the  projectile  axis  of  symmetry  in 
the  direction  of  positive  z and  considering  a particular  band  r(z)  , 
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Figure  11.  Partially  embedded  cone  with  head-on  view  of 
8 circular  band  r(z) 

one  sees  the  head-on  view  shown  in  Figure  11.  The  stress  pn  acts  on 

both  sides  of  the  band  from  iji  = 0 to  $ = ^ffiax  > and  the  dashed  arc 
represents  the  portion  of  the  band  which  is  not  embedded.  (Note  that 
the  xz  plane  is  a plane  of  symmetry  in  Figures  9-H-)  The  relation  be 
tween  i>mgx  , z , w , and  0 is  as  follows: 


z - z. 


cos  Vax  ” z tan  0 tan  to  ’ 


z^  <.  z £ z2  <_  L 


where 


x 1 + iah  0/  tan  u 


z2  ~ 1 - tan  0 tan  u 


and  z is  the  point  where  the  projectile  axis  of  symmetry  (z  axis) 

intersects  the  target  surface  (Figure  id).  The  infinitesimal  moment 

dM„_  about  the  CG  due  to  a force  dF  acting  on  any  infinitesimal 
GG 

portion  dA  of  the  cone  surface  is 


dMCG  - Y « dF 


•A.  . •• 

where  y is  the  moment  arm.  (The  moment  arm  y is  a position  vector 
from  the  CG  to  the  point  where  the  force  dF  acts.)  However,  since 

■ W 

the  xz  plane  is  a plane  of  symmetry,  the  total  vector  moment  MC(J 
about  the  CG  must  be  perpendicular  to  the  xz  plane  (i.e.,  parallel  to 
the  y axis).  Thus,  one  need  consider  only  components  of  y and  dF 
which  lie  in  the  xz  plane.  The  x and  z components  of  y are 


y (z,iji)  = z tan  0 cos 

X 


yz(z)  - z - zQ 


where 


z < L < z„ 

— U 


! m 


The  x and  z components  of  dF  acting  on  an  infinitesimal  surface 
element  dA  = z dz  difi/cos  0 are 


dF^  = -pn  tan  0 cos  (i  z dz  df 


, •i’,. 

I 


(71) 


dF  * p tan  0 s dz  d$ 
i n ~ 


Thus,  dM^  has  only  one  component  dM^  , v'nich  is  parallel  to  the 
y axis  (a- positive  moment  increases  w)  and  is  given  hy 


^CG  * yz  " 7X  * ?n  tan  #(ZG  " z 8e°2  #)z  008  ♦ dz  dlf»  (72) 


Equations  70-72  are  nov  integrated  over  the  embedded  surface  of  the  nose 
to  obtain  the  x and  z components  of  the  resisting  force,  F and 
F , and  the  total  moment  M„  : 


F = 
x 


* S' 

2p  zd  (‘ 
*n  o | 

~ tan  a>  j 


f(t)  d5 


(73) 


2 2 
F = ?rp  z tan 
z ^n  o 


2 1 a f ‘ 

’1  x J 


* (5)5  d5 

ymax  * * s 


'(71*) 


M-.„  = 


2p  z3 
rn  o 


CG  tan  w 


/ 


Sr  f(5)  d5 


**/■ 

f 


sec  0 J f (5)5  d5 

i 


(75) 


* 

Where  the  factor  "2"  arises  in  Equations  73-75  due  to  symmetry  -about 
the  xz  plane,  and 


f(?)  * (25  - 1 - x2s2) 


= (l  - ta^.2  to  tan2  0) 


5G  ~ z 


(l  + tan  0 tan  to)  * 1 


» zi  iL 


• » \ > L 


(l  - tan  0 tan  to)  * 2 


> zo  1 L 


, z2  > L 


and  \J»  , z , and  z0  are  given  "by  Equations  63-65,  The  integral 

IQ8>X  X u 

appearing  in  Equation  7^  must  he  evaluated  numerically.  The  integrals 
appearing  in  Equations  73  and  75  are  evaluated  as  follows: 


f(5)  AS 


r fz  h 

2X^  J 2X^  ^ 


X2  ^ 0 


1 . /x2e  - 1 

— arc  3in  I — 2 

x to 


, X > 0 


r 

i.:  2 

y When  X - 0 , Equations  82  and  85  reduce  to: 


i 

i 

(86)  i 


I 


where  X,Y,Z  is  the  coordinate  system  shown  in  Figures  8 and  9,  and 

F and  F are  given  by- -Equations  73  and  ?U . Likewise,  during  the 
X z 

embedding  of  the  nose,  the  rotation  of  the  projectile  about  the  CG  is 
governed  by 


Iu  = MC(,  , 0 £ zn  £ L 


where 

I = moment  of  inertia  about  CG 
a)  = angular  acceleration 

= moment  about  CG  (Equation  75) 

Ou 

The  rotational  equation  of  motion  (91 ) is  limited  to  z^  £ L , and  no 
equation  for  M^G  has  been  developed  for  z^  > L . Thus,  the  present 
theory  of  oblique  impact  is  limited  to  the  prediction  of  obliquity  w 
and  rotational  velocity  m for  0 <_  z^  £ L , i.e.,  the  theory  considers 
only  those  rotations  which  occur  during  the  embedding  of  the  nose.  A 
realistic  theory  of  deep  penetration  at  finite  angles  of  obliquity 
should  account  for  aft  body  effects  which  are  beyond  the  scope  of  the 
present  work. 


PART  III:  SOME  PARAMETRIC  IMPLICATIONS  OF  PENETRATION  THEORY 


The  Nondimensional  Equation  of  Motion 


34.  According  to  the  theory  developed  in  Part  II,  the  equation  of 
motion  for  a projectile  which  impacts  normal  to  the  surface  of  a homoge- 
neous target  is 

(M  + + W/)  <92> 

When  the  ratio  R&  is  small,  Equation  92  reduces  to 

Mii  a -’an(l>s  + W/)  • Ra<<:1  (93) 

When  the  effect  of  compressibility  is  small  (Appendix  B),  Equations  A89 
and  A93  are  approximated  by 


and 


P*-P 


Thus,  Equation  93  gives  way  to 


Mq  — -ira^  (p  + ^ pf  q^ 
n\rs  2 n1 


Equation  96  is  now  rewritten  as 


Mtj 


2 

na  p 
ors 


2 / „ . .2N 

a / 3pf  q 

-2-1  + 2_ 

%\  *• 


By  defining  the  dimensionless  quantities 


(94) 


(95) 


(96) 


(97) 


37 


.2 

u = 2_ 

U " 2 
v 

o 


(98) 


and 


= 


(99) 


where  vq  is  the  impact  velocity.  Equation  97  then  becomes 


Mq 


Mv' 


2 „ 3 dy 

Tra  p 2rra  p 

os  os 


n o 


(100) 


where 


_ . p£E2Vo  ^ 3pVo 

0=  Ps  "2ps 


(101) 


The  quantity  Gq  characterizes  the  ratio  of  maximum  inertial  resistance 
to  shear  resistance  within  the  context  of  the  cavity  expansion  and  pene- 
tration theories  presented  ir.  ohe  present  work.  Since  pg  is  primarily 
a function  of  the  target  yield  strength  Y for  the  elastic-plastic  case 
(Equation  A96),  Equation  101  suggests  another  ratio, 

Rs  H (102) 


The  ratio  Rg  is  analogous  to  the  so-called  Reynolds  number  used  in 

fluid  mechanics,  in  that  it  represents  an  order  of  magnitude  index  of 

the  ratio  of  inertial  forces  to  shear  forces  in  the  target.  Thus,  the 

ratio  R will  henceforth  be  referred  to  as  the  "solid  Reynolds  number." 
s 

35.  For  many  metals  and  earth  materials,  the  practical  range  of 
pg  is  2Y  ~ pg  £ 5Y  , so  that  at  impact  (q  = vq)  the  3range  of  Rg  is 
Gq  ~ Rg  ~ 3Gq  . Using  the  symbol  ~ to  relate  quantifies  which  are  of 
the  same  order  of  magnitude,  it  is  appropriate  to  write 


R 


s 


(103) 
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The  quantity  M$/naQps  in  Equation  100  represents  a nondimensional 
acceleration  (or  force),  integrating  Equation  100  from  y = 0 to  the 
final  nondimensional  penetration  depth'  y^  = q^/aQ  , the  following  equa- 
tion results: 


where 


yP  = sTiT  + Vo”!1  + 2?  ■ yP  ^ 25 

n a * 


u. 


\ -6R  f f ? 
1 \ a n v 

f G f 
n o/ 


1 

f G 


n o 


(10U) 


(105) 
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for  cones,  and 


(106) 


(2  - 2c  + c2)  V2c  - e2  - | (2c  - c - (1  - c)  arc  tan  - (1  - e)2  Vgc  - e2  (107) 

c2V2e  - c2 

for  ogives.  The  quantity  e is  defined  in  Part  II,  Equation  35,  and 
t,  = L/D  . After  Equation  100  is  solved  explicitly  for  u(y)  , the 
following  equations  are  obtained  for  the  acceleration  of  a projectile 
with  a conical  nose: 


= -(1  + f G u. )( 
2 n o 1 


-3R„f  (y-2t) 


a n 


, y 1 2C 


(108) 


TO.  p 

OS 


••  2 -R  f y3/4?2 

JS-‘-fc  + V.>i3«  an 

to  p 4? 

O^S 


, y i 2?  (109) 


Similar  but  more  lengthy  expressions  can  be  obtained  for  ogives.  For 
Equations  104  and  108  to  be  applicable,  it  is  obvious  that  the  projec- 
tile must  penetrate  at  least  one  nose  length  L . After  having  inte- 
grated the  approximate  equation  of  motion  (100)  in  nondimensional  form, 
it  is  now  possible  to  examine  explicitly  the  roles  of  the  parameters 
G , R , R . and  £ = L/D  . Furthermore,  while  the  added  mass  term 

U 5 cl 
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ira  p.B.  has  been  neglected  in  Equations  93-109,  the  ratio  R still 

O X*  X & 

appears  as  a parameter  in  the  nondimensional  final  depth  of  penetration 
y - 9p/aQ  , Equation  lOU. 

The  Shape  of  the  Deceleration  Curve 

3 6.  Probably  the  most  stringent  test  of  the  present  penetration 
theory  lies  in  the  comparison  of  predicted  deceleration  versus  depth  or 
time  with  experimental  observations  and  with  the  results  of  detailed 
two-dimensional  finite-difference  calculations,  where  available.  For 
the  moment,  however,  attention  is  focused  on  the  shape  of  the  decelera- 
tion curve  in  order  to  make  a qualitative  comparison  with  experimental 
results.  Projectiles  with  conical  noses  are  amenable  to  such  an  exam- 
ination because  of  the  simplicity  of  Equations  108  and  109,  but  it  is 
emphasized  that  similar  results  can  be  obtained  for  ogives. 

37-  The  deceleration  curve  has  the  same  general  shape  whether 
depth  or  time  is  chosen  as  the  abscissa.  The  overall  shape  of  the  curve 

is  strongly  dependent  on  G ~ R and  on  X,  ~ L/D  . This  is  illus- 

o s 

trated  in  Figures  12  and  13.  The  increasing  portion  of  the  curve  cor- 
responds to  the  embedding  of  the  nose,  and  the  decreasing  portion  cor- 
responds to  a fully  embedded  nose.  Examination  of  experimental 
resultsl!,:,1'i,16  reveals  that  many  observed  deceleration  curves  tend  to 
be  somewhat  flat  after  the  nose  becomes  embedded  in  the  target,  thus 
having  the  general  shape  of  the  curves  in  Figure  12  (at  least  for  earth 
penetration),  but  the  rapid  decay  shown  in  Figure  13  is  rarely  seen. 

The  present  penetration  theory  leads  to  reasonably  accurate  predictions 

of  final  penetration  depth  for  0 < R <,  100  , but  only  the  average 

s 

value  of  the  deceleration  is  correctly  predicted  when  R is  signifi- 

s 

cantly  greater  than  one.  Moreover,  the  theory  apparently  overpredicts 
the  instantaneous  target  inertial  resistance  during  the  early  part  of 
the  penetration  process,  and  the  magnitude  of  this  overprediction  in- 
creases as  R increases  from  one.  The  net  result  is  that  the  theory 
s 

tends  to  underpredict  final  penetration  depths  when  R is  very  large 

s 

at  impact. 


The  Relative  Effects  of  Nose  Shape  and  Frontal  Loading 
on  Final  Depth  of  Penetration 


38.  In  the  context  of  the  present  penetration  theory,  the  effect 
of  nose  shape  is  coupled  with  the  effects  of  R&  and  Gq  ~ Rg  , so  that 
the  final  depth  of  penetration  will  not-  he  proportional  to  some  function 
of  nose  shape  alone.  The  nose  performance  ratio  Rr  is  now  defined  as 

R = — (110) 

n Sb 

where  q^  is  the  final  depth  of  penetration  for  a projectile  with  a 
given  nose  shape,  and  q^  is  the  corresponding  quantity  for  an  other- 
wise identical  blunt  cylinder  {?  = 0).  Combining  Equations  10U  and  110 
leads  to 


In  (1  + f Gu  ) + 6cf  R 
o _ n o 1 n a 

n = f In  (1  + G ) 
n o 


(111) 


where  u^  is  given  by  Equation  105.  In  Figures  14-17,  Rn  is  plotted 
versus  £ for  cones  and  ogives  with  R&  and  Gq  as  parameters.  Fig- 
ures l4  and  16  correspond  to  shallow  penetration  depths  on  the  order 
of  a few  nose  lengths,  and  Figures  15  and  17  correspond  to  depths  sev- 
eral times  deeper  than  Figures  l4  and  16.  The  values  R =0.1  and 

B. 

R = 0.001  were  chosen  to  conform  approximately  to  the  upper  and 

€L 

lower  limits  of  R in  penetration  tests  conducted  by  Sandia  Labora- 
15  16  ^ 

tories,  ’ and  the  value  R = 0.01  represents  an  intermediate  value. 

a 2 

For  comparison  Young's  empirical  correlation  is  included  in  Fig- 
ures 14-17.  Young's  linear  relation  between  ? and  final  depth  of 
penetration  was  deduced  from  earth  penetration  data  for  a variety  of 

projectiles  and  impact  conditions  and  represents  something  of  an  overall 

2 

average  nose  shape  effect. 

39.  The  present  theory  results  in  a nose  shape  effect  which  be- 
comes more  pronounced  as  R increases.  This  simply  means  that,  for 

£L 

a given  projectile  diameter,  a sharp  nose  is  less  effective  for  a heavy 
projectile  than  for  a light  projectile  as  far  as  depth  of  penetration  is 


42 


Wf/r— 


legend 


matwwaK:  EQUATION  III, 
' R.  * o.t  . a,  = 3.0 


R*  = 0.1  , G0  = 5.0 
R*  = 0*0)  . G_  = 0.1 


R.  = 0.01  . Gfl  = 0J 

* R,  - 0.001.  G#  = 0.03 

6 R,  = 0.001,  G.  = 0.05 

QmEJ  £MP!R'CAL  correlation. 


igure  16.  Nose  performance 


ratio  for  ogives,  "shallow" 


penetration 


LEGEND 


SSSLhS  mesent  work.  equ«t,on 

' R.  * »•'  . c.  . io 

Ra  “ °*)  • G.  = too 


« . - 10 

Ra  = °*)  . G0  = ioo 

Ra  = 0.01  , G„  = i 

Ra  = 0.01  . G„  = ,o 


c » “w*  • t.0  = 10 

r Ra  = 0-001,  Ge  = , 

Ra  = 0.001,  G#  = ,0 

Soever  empirical  correlation 
young2 


l/d 


Figure  17.  Ifose  performance 


ratit.  for  ogives,  "deep"  penetration 


concerned.  The  present  theory  also  infers  that  nose  shape  is  more 
effective  in  shallow  penetration  where  the  embedding  of  the  nose  con- 
stitutes a significant  portion  of  the  penetration  process. 

UO.  There  have  been  various  attempts  to  ascertain  the  relation 

p 

between  final  penetration  depth  and  frontal  loading  W/tra^  = W/A  , where 
W is  the  projectile  weight  and  A is  the  projected  frontal  area. 
Young's  empirical  correlation,  for  example,  infers  that  final  penetra- 
tion depth  is  proportional  to  the  square  root  of  W/A  . The  present 
theory,  on  the  other  hand,  yields  an  approximate  linear  relation  between 
final  (dimensional)  penetration  depth  and  W/A  for  depths  greater  than 
one  nose  length  (Equation  104) . 

Parametric  Differences  Between  Theoretical  and 
Empirical  Results  for  Homogeneous  Targets 

kl.  Within  the  context  of  the  present  penetration  theory,  it  is 
evident  that,  for  the  penetration  of  homogeneous  targets,  final  pene- 
tration depth  is  approximately  proportional  to  W/A  and  to  the  natural 
log  of  the  square  of  the  impact  velocity.  Furthermore,  the  effect  of 
nose  shape  is  coupled  with  the  effects  of  W/A  , impact  velocity,  and 

target  constitutive  properties.  Young's  empirical  equation  fo”  final 
2 

penetration  depth,  however,  indicates  that  the  final  depth  is  propor- 
tional to  nose  length,  impact  velocity,  and  the  square  root  of  W/A 
regardless  of  target  properties.  Moreover,  the  relative  penetrability 
of  the  target  is  characterised  by  a single  multiplicative  constant  which 
must  be  derived  from  penetration  experiments.  Although  Young's  equation 
may  correctly  reflect  the  relative  effects  of  impact  velocity,  nose 
shape,  and  W/A  , it  represents  a statistical  correlation  of  experi- 
mental observations  and  offers  little  explanation  in  terms  of  underlying 
physical  phenomena.  On  the  other  hand,  the  present  penetration  theory 
is  based  on  certain  assumptions  about  the  underlying  physical  processes 
but  produces  trends  somewhat  different  from  Young's  equation.  Even 
though  the  theory  is  "rationally”  based,  it  may  not  be  sufficiently 


realistic;  a more  realistic  theory  might  result  in  trends  similar  to 
those  exhibited  “by  Young's  equation. 

Penetration  of  Layered  Targets 

k2.  The  relatively  complex  treatment  of  the  layered  target  in 
Part  II  is  important  for  predicting  projectile  loading  when  the  initial 
layer  thickness  is  comparable  to  the  projectile  hose  dimensions  and/or 
•when  the  layers  have  extremely  different  properties.  On  the  other  hand, 
when  the  layers  axe  thick  and  have  fairly  similar  properties,  final 
depth  predictions,  which  approximate  the  results  of  the  layered  target 
theory  can  "be  obtained  simply  by  using  the  theory  for  homogeneous  tar- 
gets and  changing  the  material  property  input  when  the  projectile  reaches 
a layer  interface. 
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PART  IV:  APPLICATIONS  OF  PENETRATION  THEORY 
Introduction 

* 43.  The  penetration  theory  presented  in  Part  II  is  relatively 
simple  compared  with  two-dimensional  finite-difference  analyses.  In 
order  to  demonstrate  in  part  the  applicability  of  the  theory  as  an  engi- 
neering tool,  theoretical  results  will  he  compared  with  two-dimensional 
finite-difference  solutions,  empirical  results,  and  experimental  obser- 
vations. However,  before  these  comparisons  are  made,  some  comments  are 
in  order: 

a.  Although  the  ratio  R « 1 in  all  the  examples  to  be 

a 3 

presented,  the  added  mass  term  Tta^p^B^  will  be  retained 
as  an  upper  limit  on  the  effect  of  target  acceleration. 

b.  It  will  be  shown  that  the  solid  Reynolds  number  R 

' s 

provides  a simple  rule  of  thumb  which  indicates  whether 
or  not  the  present  theory  is  appropriate  for  a particular 
penetration  problem. 

c.  Material  properties  and  impact  conditions  are  given  in 
both  U.  S.  customary  and  metric  (Si)  units  in  the  text, 
but  graphic  presentation  of  results  is  given  in  U.  S. 
customary  units  to  expedite  comparison  with  other  work. 

Choice  of  Material  Properties 

44.  The  cavity  expansion  theory  presented  in  Appendix  A and  the 
penetration  theory  developed  in  Part  II  are  restricted  to  materials 
which  exhibit  bilinear  behavior  in  shear,  as  shown  in  Figure  A2.  More- 
over,. the  assumption  of  rate-independent  shear  behavior  is  implicit  in 
the  theoretical  development.  Real  materials  exhibit  neither  bilinear 
shear  behavior  nor  rate-independence,  and  the  guidelines  to  be  followed 
when  the  present  penetration  theory  is  applied  to  real  targets  are: 

(a)  bilinear  stress-strain  curves  should  be  chosen  so  as  to  approximate 
the  total  area  under  experimentally  obtained  stress-strain  curves,  up  to 


the  point  where  complete  failure  occurs  (i.e. , where  the  slope  becomes 
negative);  and  (b)  when  predictions  are  to  be  made  prior  to  the  penetra- 
tion of  an  untried  target,  it  is  best  to  make  two  sets  of  calculations 
corresponding  to  the.  upper  and  lower  bounds  on  Y , E , and  . Lower 
bound  properties  can'  be  obtained  from  the  results  of  static  tests,  while 
upper  bound  properties  must  be  estimated  from  the  results  of  dynamic 
tests. 


Fenetration  of  Rock 


U5.  Thigpen  has  performed  two-dimensional  finite-difference  cal- 
culations analyzing  the  penetration  of  a nonrigid  projectile  into  Madera 
limestone  and  welded  tuff.  The  present  penetration  theory  is  applied  to 
the  same  problems  for  a rigid  but  otherwise  identical  projectile.  Fol- 
lowing Thigpen,  the  targets  are  idealized  as  compressible  elastic- 
plastic  continua  of  the  von  Mises  type  (compressive  yield  strength  Y 
independent  of  pressure).  For  the  range  of  dynamic  pressures  considered 
(P  < Pvq),  the  relation  between  pressure  ana  density  is 


where  P 

ave 

dix  B) , and 
constant  for 


+ <m> 

is  average  dynamic  pressure  in  the  plastic  zone  (Appen- 
CQ  is  initial  sound  speed.  The  modulus  E is  essentially 
the  range  of  pressures  examined.  For  the  limestone, 


p = 168  pcf*  (2.69  gm/cra^) 

E ss  3.15  x 106  psi  (217  kbar) 

Y = 1.37  * 101*  psi  (0.945  kbar) 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  is  presented  on  page  4. 


vatsameot** 


and 


For  the  tuff. 


* Ji  o 

C = 1'.12  x 10  fps  (3.4l  x 10J  m/sec) 


p = 115  pcf  (1.84  gm/cnr*) 

E — 1.09  x 106  psi  (75. 2 kbar) 


and 


Y - 5-51  x 105  psi  (0.380  khar) 


C = 6.72  x io3  fps  (2.05  x 103  m/sec) 


For  limestone  penetration, 


1 


and 


W = 674  lb  (306  kg) 
D = 8 in.  (20.3  cm) 
L/D  = 3.0  (9.25  CRH) 


vq  = 570  fps  (174  m/sec) 


| 

I 


For  tuff  penetration, 

W = 1000  lb  (454  kg) 

D = 9 in.  (22.9  cm) 

L/D  = 2.4  (6.0  CRH) 
and 

vo  = 695  fps  (212  m/sec) 

46.  The  predictions  of  the  present  penetration  theory  are  com- 
pared with  Thigpen's  results  in  Figures  10  and  19 . Thigpen's  decelera- 
tion curves  correspond  to  the  motion  of  a particle  on  the  projectile 
axis  of  symmetry;  in  Figure  18  the  particle  is  located  2 ft  (6l  cm)  aft 
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Figure  l8.  Motion  curves  for  penetration  into  limestone 


Figure  19.  Motion  curves  for  penetration  into  tuff 
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of  the  nose  tip,  and  in  Figure  19  the  particle  is  located  at  the  "base  of 
the  nose.  The  oscillation  in  Thigpen's  deceleration  curves  is  due  to 

7 

nonrigidity  of  the  projectile.  It  is  important  to  observe  that  in 

these  calculations  the  target  shear  resistance  is  very  significant  since 

the  values  of  Rg  at  impact  are  0.86  and  2.17  for  the  limestone  and 

tuff,  respectively.  Moreover,  the  agreement  with  Thigpen's  results 

apparently  provides  justification  for  using  p to  characterize  average 

s 

shear  resistance  to  penetration,  at  least  for  homogeneous  targets  and 

low  values  of  R . 

s 


Penetration  of  Concrete 


47.  Another  test  of  the  present  penetration  theory  lies  in  the 
comparison  of  theoretical  predictions  with  those  of  empirical  penetra- 
tion equations.  A fairly  accurate  empirical  equation  for  the  penetra- 
tion of  armor  piercing  (AP)  and  semi-armor  piercing  (SAP)  bombs  into 

17 

reinforced  concrete  is  given  by 


°*P 


222WD0,21^ 

aVy 


+ 


0.5D  + 15% 


(113) 


with  q^  in  inches,  W in  pounds,  A in  square  inches,  D in  inches, 

Y in  psi,  and  v in  fps.  Equation  113  is  now  used  to  predict  the 

penetration  of  a 2000-lb  (907 -kg)  AP  bomb  into  reinforced  concrete, 

6 

with  Y = 5000  psi  (0.345  kbar),  E = 4.10.  x 10  psi  (283  kbar),  and 

3 

p = 150  pcf  (2.4  gm/cm  ) . For  the  penetration'  theory  predictions,  the 

modulus  E.  is  taken  to  be  zero,  and  the  value  of  YNvis  the  static 

* X 

value  (5000  psi)  as  in  Equation  113.  In  addition,  the  prbasure-density 

curve  in  Figure  20  is  used  according  to  the  procedure  outlinecKin  Appen- 
dix B and  represents  typical  behavior  of  concrete  under  pressure. '\The 
bomb  has  D = 13.5  in.  (34.3  cm)  and  L/D  = 1.12  (1.5  CRH)  . Normal^ 
impact,  velocities  range  from  500  to  1500  fps  (152  to  457  m/sec).  Two 
sets  of  theoretical  predictions  are  made,  corresponding  to  the  compress- 
ible and  incompressible  cases.  These  are  compared  with  the  results 
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Figure  20.  Typical  pressure-density  relation  for 
concrete  (l  ksi  = 1000  psi) 

of  Equation  113  in  Figure  21.  Since  the  penetration  theory  predictions 
are  made  us^ng  the  static  value  of  Y , they  should  he  regarded  as  upper 
bounds  on  predicted  penetration  depth.  The  results  in  Figure  21  corre- 
spond to  impact  values  of  R in  the  range  1 <,  R < 15  , in  which  the 
inertial  resistance  of  the  target  begins  to  predominate  over  the  shear 
resistance.  Moreover,  the  agreement  with  Equation  113  provides  some 

justification  for  using  the  theory  in  this  range  of  values  of  R . 

s 
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EMPIRICAL  RESULTS,  EQUATION  1 13 

PRESENT  WORK-COMPRESSIBLE 
PRESENT  WORK-INCOMPRESSIBLE 


IMPACT  VELOCITY,  FPS 


Figure  21.  Penetration  depth  versus  impact  velocity 
for  penetration  of  concrete  by  a 2000-lb  AP  bomb 


Penetration  of  Soft  Earth 


1*8.  As  an  example  of  penetration  into  soft  earth,  the  normal  im- 
pact of  a terminal  delivery  vehicle  (TDV)  into  moist,  sandy  clay  is  con 
sidered.  The  target  is  one  of  the  range  areas  at  Jefferson  Proving 


I 


Ground,  Indiana.1^  The  TDV  has  a total  weight  of  17.3  lb  (7*85  kg),  “but 
when  the  full  length  of  the  projectile  is  embedded  in  the  target,  the 
projectile  base  separates,  remaining  at  the  target  surface  and  reducing 
the  projectile  weight  to  13.6  lb  (6.17  kg).  The  diameter  of  the  pro- 
jectile is  3 in.  (7.62  cm),  the  total  projectile  length  is  17.7  in. 

(45.0  cm),  and  the  nose  is  an  ogive,  L/D  = 1.45  (2.35  CRH).  For  the 

target,  Y = 50  psi  (3.45  bar),  E = 3000  psi  (207  bar),  and  E.  = 0, 

19 

obtained  from  dynamic  unconfined  compression  tests.  The  undisturbed 

O O 

target  density  p = 3.88  slugs/ft  (2.00  gm/cm  ).  Experimental  impact 
velocities  range  from  350  to  610  fps  (107  to  186  m/sec)  corresponding  to 
65  ~ ~ 200  . Since  no  pressure-density  curve  is  available  for  the 

target,  two  sets  of  calculations  are  presented,  corresponding  respec- 
tively to  the  incompressible  case  and  the  compressible  case  for 

= l.lp  (probably  representing  a high  upper  limit  on  the  effect  of 
compressibility).  The  theoretical  predictions  are  compared  with  experi- 
mental results  in  Figure  22.  Since  the  experimental  data  are  quite 
scattered,  it  is  difficult  to  come  to  any  definite  conclusions  from  this 
comparison.  However,  within  the  framework  of  the  theory,  the  incom- 
pressibile  calculation  is  probably  the  more  realistic  of  the  two, 
although  the  predicted  depths  are  somewhat  low  compared  to  the  experi- 
mental results.  Viewed  in  this  light,  the  results  seem  to  indicate  that 

R ~ 100  may  represent  an  approximate  upper  limit  on  the  range  of  ap- 
s 

placability  of  the  present  penetration  theory. 

< 

' Perforation  of  a Metal  Slab 

49.  Perhaps  the  most  stringent  test  which  can  be  applied  to  the 

* v 

penetration .theory  for  layered  targets  developed  in  Part  II  is  a situa- 
tion in  which  the  first  layer  is  thin  with  extremely  different  proper- 
ties from  the  second  layer.  Such  a situation  exists  when  a steel  pro- 
jectile with  a diameter  of  0.3  in.  ( 0.762  cm)  strikes  a 606I-T6 
aluminum  slab  with  a thickness  of  1 in.  (2.54  cm).  This,  of  course, 

represents  a case  where  the  first  layer  is  aluminum  and  the  second  layer 

20 

is  air.  Using  the  HEMP  code,  Wilkins  has  performed  a two-dimensional 
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Figure  22.  Penetration  depth  versus  impact  velocity  for  pene- 
tration of  soft  earth  by  a terminal  delivery  vehicle 

finite-difference  analysis  of  the  perforation  of  1-in. -thick  aluminum 

slab  by  a 0.3-in.-diam  steel  projectile.  The  projectile  has  a conical 

nose,  L/D  = 1 , a weight  of  0.0183  lb  (8.30  gm),  and  a normal  impact 

velocity  of  2758  fps  (840  m/sec).  For  the  aluminum,  p = 5.25  slugs/ft3 

(2.7  gm/crn3),  Y = 43,500  psi  (3  kbar),  S = 1.03  * 107  psi  (710  kbar), 

and  E = 0 . The  pressure-density  relation  is  Equation  112  with 
2 7 

pCQ  = 1.24  x 10'  psi  (775  kbar).  The  properties  of  air  are  considered 
negligible  compared  to  those  of  aluminum.  According  to  Wilkins,  the 
projectile  perforates  the  slab  with  an  exit  velocity  equal  to  8l  percent 
of  the  impact  velocity.  This  represents  a projectile  momentum  loss  of 
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19  percent  and  a kinetic  energy  loss  of  3^.^  percent.  The  present  pene- 
tration theory  for  layered  targets  has  been  applied  to  the  same  problem, 
arid  the  results  are  presented  in  Figures  23-26.  The  crucial  test  of  the 
theory  lies  in  its  ability  to  predict  projectile  kinetic  energy  loss 
during  the  perforation  process i The  incompressible  theory  yields  an 
energy  loss  of  percent,  and  the  compressible  theory  reduces  this  to 
52  percent,  essentially  the  ssme  answer.  Thus,  the  present  theory  over- 
predicts Wilkins*  value  for  the  energy  loss  by  a minimum  factor  of  1.51. 
Two  possible  sources  of  error  are  immediately  apparent:  First,  as  is 

pointed  out  in  Part  III,  paragraph  37,  the  theory  apparently  overpre- 
dicts the  initial  target  inertial  resistance  when  the  value  of  R is 

s 

significantly  greater  than  one  at  impact  (in  this  case,  R =6.29  at 

s 

impact).  Second,  the  target  is  always  assumed  to  be  in  a compressive 
state  of  stress;  however,  it  is  well  known  that  because  the  backside  of 
the  slab  is  a free  surface,  a finite  region  in  the  slab  may  be  in  a ten- 
sile state  of  stress  during  the  perforatioz^process,  thus  reducing  the 
overall  penetration  resistance. 

50.  Within  the  framework  of  the  theory  as  it  stands,  there  is 
little  that  one  pan  do  about  the  suspected  region  of  tensile  stress,  and  ^ 
this  may  represent  a permanent  limitation  on  the  theory.  On  the  other 
hand,  the  question  of  inertia  is  a bit  more  straightforward.  Accord- 
ingly, a third  calculation  is  presented  in  which  the  target  is  assumed 
incompressible,  and  the  inertial  terms  in  the  projectile  equation  of 
motion  are  discarded  altogether.  For  this  case,  the  predicted  energy 
loss  is  32  percent,  slightly  underpredicting  Wilkins*  value  of  3^.1*  per- 
cent.  This  does  not  necessarily  mean  that  all  the  target  resistance  is 
due  to  shear.  However,  it  may  indicate  that  the  theory  does  a better 
job  of  predicting  target  shear  resistance  than  predicting  target  inertial 
resistance.  In  any  case,  it  is  emphasized  that  the  present  theory  is 
developed  not  as  a perforation  theory,  per  se,  but  as  a penetration 
theory  for  layered  targets.  For  most  cases  of  practical  interest,  the 
layers  will  have  comparable  densities  (i.e.  comparable  inertial  resis- 
tance) even  though  other  properties  may  bt.  radically  different.  Since 
the  application  of  the  layered-target  theory  to  a perforation  problem 
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Figure  24.  Projectile  kinetic  energy  versus  depth  for 
perforation  of  1-in. -thick  aluminum  slab 
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Figure  25.  Projectile  velocity  versus  time  for  perforation 
of  1-in. -thick  aluminum  slab 


represents  a very  extreme  case  and  yields  fair  results,  the  theory  prob- 
ably represents  a reasonable  extension  of  the  penetration  theory  for 
homogeneous  targets. 


Oblique  Impact  Against  Hard  Earth 

51.  At  the  present  time,  no  data  for  oblique  impact  and  ricochet 
have  been  found  which  are  appropriate  for  comparison  with  the  oblique 
impact  theory  developed  in  Part  II,  i.e. in  the  experimental  results 
found  to  date,  there  is  not  enough  information  given  about  projectiles 
(e.g. , CG  location  and  moment  of  inertia)  to  effect  a reliable  test  of 
the  theory.  Nevertheless,  an  application  of  the  theory  to  the  oblique 
impact  of  a TDV  against  a hard  earth  target  is  made  without  experimental 
verification  and  is  presented  as  a sample  calculation.  The  projectile 
characteristics  are 


W = 15  lb  (6.80  kg) 

I = 0.0975  slug-ft2  (0.132  kg-m2) 


The  target  is  equivalent  to  frozen  sand  found  at  Camp  McCoy,  Wisconsin, ^ 
with  . ' 

p = = 3.21  slugs/ft^  (1.66  gm/cm^) 


X=  958  psi  (66.1  bar) 


E = 6.25  * lCT  psi  (It. 31  kbar) 


Et  = ° 

(These  properties  have  teen  previously  used  in  the  analysis  of  TDV 
impact  and  penetration  at  Camp  McCoy. The  initial  angle  of  obliguity 
is  1*5  deg  (0,76  rad),  and  the  impact  velocities  considered  are  300  fps 
(91.1*  m/sec)  and  600  fps  {183  m/sec).  Initial  and  final  projectile 
positions  are  illustrated-  for  each  impact  velocity  in  Figures  27  and  28, 
demonstrating  the  dependence  of  vehicle  rotation  bn  impact  velocity.  In 
the  300-fps  case,  the  projectile  undergoes  a rotation  of  15.1  deg,  and 
the  downward  motion  of  the  CG  is  stopped  before  the  nose  is  fully  em- 
bedded. In  the  600-fps  case,  the  projectile  undergoes  a rotation  of 
li.17  deg  during  the  embedding  of  the  nose.  The  curve  connecting  the  CG 
locations  is  the  CG  trajectory. 
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PART  V;  CONCLUSIONS  AND  RECOMMENDATIONS 

52.  The  results  presented  in  Part  IV  stand  as  partial  evidence 

that  the  cavity  expansion  theory  can  "be  used  as  an  approximate  basis  for 

the  development  of  a projectile  penetration  theory,  particularly  when 

the  effect  of  shear  is  predominant  (i.e.,  for  low  values  of  the  solid 

Reynolds  number  R ).  However,  it  is  emphasized  that  the  analogy  he- 
s 

tween  spherical  cavity  expansion  and  projectile  penetration  is  somewhat 
fortuitous  ,■  and  the  present  penetration  theory  probably'  reflects  only 
the  gross  features  of  the  projectile-target  interaction.  Thus,  the 
acquisition  of  insight  into  the  mechanics  of  penetration  demands  the 
application  of  whatever  means  are  available,  including  empirical  methods 
and  two-dimens ional  finite-difference  analyses. 

53.  For  estimates  of  penetration  depth,  the  applicability  of  the 

present  penetration  theory  appears  to  lie  in  the  approximate  range 

0 < R $.  100  ; but  the  upper  bound  is  offered,  as  a conservative  rule  of 
s 

thumb  and  not  as  an  absolute  limit.  (Results  of  experiments  at  high 

values  of  R are  needed  in  order  to  establish  a more  realistic  range . ) 
s 

However,  it  is  important  to  note  that  as  R ->  0 the  problem  of  dynamic 

s 

penetration  reduces  to  that  of  a static  punch,  in  which  friction  may  be 
important.  On  the  other  hand,  as  03  shear  resistance  becomes 

insignificant  such  that  the  projectile-target  interaction  may  reduce 
essentially  to  a fluid  dynamics  problem. 

54.  Although  some  footholds  have  been  established,  the  develop- 
ment of  the  penetration  theory  is  still  in  an  embryonic  stage;  some  sug- 
gestions for  additional  work  are  listed  below: 

a.  At  the  present  time,  the  role  of  target  inertia  in  the 
projectile-target  interaction  is  only  partially  under- 
stood, and  the  degree  of  contact  between  the  projectile 
nose  and  the  target  has  not  been  established.  Conse- 
quently, further  investigation  is  necessary  to  establish 
the  distribution  of  dynamic  pressure  on  the  projectile- 
target  contact  surface  as  a function  of  projectile  veloc- 
ity, nose  shape,  and  target  properties. 
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Although  the  cavity  expansion  shear  resistance  term  pg 
seems  to  represent  a reasonable  first  approximation  for 
average  normal  stress  at  the  projectile  surface  due  to 
'Shear  in  the  target,  some  effort  should  be  directed 
toward  establishing  normal  and  tangential  stress  compo- 
nents on  the  projectile  surface  due  to  target  shear 
behavior.  Furthermore,  the  variation  of  these  quantities 
with  nose  shape,  penetration  depth,  and  confining  pres- 
sure should  be  determined. 

The  present  oblique  impact  theory  is  limited  to  the 
embedding  of  the  projectile  nose,  with  all  rotation 
occurring  in  a single  plane.  Moreover,  this  work  has  not 
been  verified  experimentally.  Future  efforts  should 
first  assess  the  accuracy  of  the  existing  theory.  Then, 
if  necessary,  the  theory  can  be  extended  to  include  aft 
body  effects  and  motion  in  three  dimensions.  Eventually, 
it  may  be  possible  to  formulate  a workable  theory  for 
deep  penetration  at  finite  angles  of  attack  and  obliquity 
However,  it  is  emphasized  that  overcomplication  is  as 
undesirable  as  oversimplification,  and  modification  of 
the  theory  should  proceed  one  step  at  a time  to  avoid  the 
incorporation  of  superfluous  details . 
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APPENDIX  A:  SPHERICAL  CAVITY  EXPANSION  IN 

A CONCENTRICALLY  LAYERED  MEDIUM 


Introduction 


1.  This  appendix  presents  an  extension  of  cavity  expansion  the- 

'8*  10 

ory  developed  by  Bishop,  Hill,  and  Mott,  Goodier,  and  Hanagud  and- 
12 

Ross.  Goodier* s treatment  dealt  with  the  dynamic  expansion  of  a 
spherical  cavity  in  an  infinite,  incompressible,  elastic-plastic, 
strain-hardening  medium.  By  means  of  a locking  approximation  in  which 
the  plastic  zone  is  assumed  to  be  uniformly  compressed,  Hanagud  and  Ross 
modified  the  theory  to  account  for  the  case  when  the  effects  of  compres- 
sion predominate  over  the  effects  of  shear.  However,  the  inertial  terms 
in  the  Hanagud-Ross  treatment  do  not  reduce  to  those  of  Goodier  for  the 
incompressible  case.  The  present  work  eliminates  this  contradiction  and 
also  accounts  for  the  effect  of  concentric  layering. 


Problem  Formulation 


2.  In  Figure  Al,  a spherical  cavity  of  radius  a(t)  is  surrounded 


LEGEND 

a = CAVITY  RADIUS 
b = RADIUS  OF  PLASTIC  FRONT 
h = INTERFACE  RADIUS 
t = TIME 


LAYER 
" NO.  1 


LAYER 
NO.  2 


Figure  Al.  Spherical  cavity  surrounded  by  concentric  layers 


* Raised  numbers  refer  to  similarly  numbered  items  in  "References"  on 
pp.  65-66  at  end  of  main  text. 
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by  two  concentric  layers  of  different  materials.  Both  materials  are 

9-12 

elastic-plastic  with  linear  strain-hardening,  as  illustrated  in 
Figure  A2. 


Figure  A2.  Material  behavior  in  shear 

Both  materials  are  incompressible  in  the  elastic  and  plastic  states,  but 

transition  from  the  elastic  to  the  plastic  state  is  accompanied  by  a 

finite  volumetric  strain  e.  . This  quantity  is  known  as  a locking 
11  12  * 

strain  * and  is  introduced  in  order  to  obtain  a first  approximation 

to  the  effect  of  compressibility.  Although  is  taken  to  be  constant 

in  the  present  analysis,  it  must  be  a function  of  pressure  to  have  any 

physical  meaning.  However,  the  variation  of  with  pressure  does  not 

alter  the  analysis  or  conclusions  of  this  appendix,  and  a discussion  of 

the  locking  epproximation  is  given  in  Appendix  B. 

3.  The  elastic  and  plastic  zones  are  separated  by  a spherical 

plastic  front  of  radius  b(t)  (Figure  Al).  Since  b(t)  is  the  only 

radial  location  where  compression  occurs,  it  also  represents  a shock 

front  where  density  and  stress  are  discontinuous.  The  equation  of  mo- 
ll 12 

tion  for  spherical  symmetry  5 at  r 5*  h on  either  side  of  b(t)  is 
then 


A2 


3<J  „ 

»r-  + l 

3r  t 


(dr  - °s>  * <‘(k  + vfe)v  5 p it  ’ b<r< 


r 2 

IT  + F (or  - ff8>  ' p* 


/i-  + v s_\ 

^3t  VA  3r  y 


V£  = P£  Dt  * 


a < r < b 


where 

ct  = radial  stress 
r 

Og  = tangential  stress 
r = radial  coordinate 
0 = polar  coordinate 

v = radia'l  velocity  in  the  elastic  region 
v^  = radial  velocity  in  the  locked  plastic  region 
U.  The  material  is  required  to  he  stress-free  as  r 00  , and  the 
discontinuity  in  at  r = b is  defined  by 


. „<+)  „(-) 

4ob  a «b  - % 


where  o'  is  the  limiting  value  of  a as  r ■+  b in  the  elastic 
b (-)  r 

region,  and  cr^  is  the  corresponding  quantity  in  the  locked  plastic 
region.  The  equation  of  motion  is  then  integrated  with  respect  to  r 
from  r = 00  to  r = a , accounting  for  the  discontinuity  at  r = b , so 
that  the  compressive  normal  stress  at  the  cavity  surface  is  given  by 


p(t)  = -or(a,t)  = + Im(a,b)  + Im(b,“)  + Ig(a,»)  (A4) 


where  the  momentum  integrals  Im(a,b)  and  I^b,”)  are  given  by 


(a,b)  = J 


p»Drdr  a”4  yb 


and  the  general  shear  integral  is  defined  by 


(A6) 


Is^rl,r2^  = 2 


Inertial  Terms 


5.  Attention  is  now  focused  on  the  terms  A a,  , I (a,b)  , and 

b m 

Im(b,»)  which  embody  the  inertial  effects  in  Equation  A4.  Consider 
first  the  stress  discontinuity  Aa^  at  r = b(t)  . Conservation  of 
■mass  at  b(t)  requires  that^"’^11" 


P£(b  - vz)  = p(b  - v)  at  r = b 


(A7) 


where  a dot  denotes  time  derivative.  Equation  A7  leads  to 


v - ab 

v = ttit  at  r - b 


!A8) 


where  a is  the  locking  compressibility  defined  by 


a = 1 - = l-  e 

Po 


(A9) 


and 


e.  =ln  ^ 

1 W/ 


Conservation  of  momentum  at  b(t)  requires  that 


(A10) 


Aafe  = p£v  (b  - v^)  - pv(b  - v)  at  r = b 


(All) 


Incorporating  Equation  A8  in  Equation  All  yields 


A4 


“W  - Y1  , 

Agbs,"l-d  afc  r = * 


(A12) 
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Conservation  of  mass  in  the  locked  plastic  region  gives  rise  to  ’ 


f,(t) 

Vr't)  " ~T~ 

r 


(A13) 


t At) 

v£(a,t)  * A = — ^ — at  r = a(t) 
aT 


(Aik) 


Thus, 


= a a 


Likewise,  conservation  of  mass  in  the  elastic  region  demands  that 


(A15) 


v(r,t)  = ^1^- 

r 


(A16) 


Substituting  Equations  A13,  A15,  and  Al6  into  Equation  A 9 and  evaluating 
at  r = b(t)  obtains 


f(t)  = 


2.  ,2/ 

a a - abl? 

1 - a 


(A17) 


Consolidating  the  definition  of  D/Dt  given  in  Equations  A1  and  A2  with 
Equations  A13  and  Al6  results  in 


Dv  f 2f  . £ 

l _ _l  Dv,  f_  2f 

Dt  " 2 ~ 5 ana  Dt  " 2 " 5 

r r r r 


(A18) 


By  combining  Equation  Al8  with  Equation  A5,  the  momentum  integrals 
become 


A5 


6.  So  far,  the  discussion  has  avoided  explicit  consideration  of 

the  layer  interface  h(t)  . However,  the  location  of  b(t)  relative  to 

h(t)  determines  the  character  of  the  momentum  integrals  I (a,b)  and 

m 

Im(b,“)  as  well  as  the  stress  discontinuity  A cr^  . Denoting  first  and 
second  layer  properties  by  the  subscripts  1 and  2,  respectively,  examine 
the  situation  when  b < h in  Figure  A3  below. 

LEGEND 

a = CAVITY  RADIUS 
b = RADIUS  OF  PLASTIC  FRONT 
h = INTERFACE  RADIUS 

Figure  A3.  Relative  location  of  cavity  surface,  plastic 
front,  and  layer  interface  when  b < h 

For  this  case, 


I (a,b)  = p 


b < h 


(A21) 


- P1  / (^2  - 2 ®r  + P2  jf  - 2 *5)  •’ 


b < h (A22) 


Evaluating  integral  in  Equations  A21  and  A22  and  employing  Equation  A13 
in  Equation  A12  result  in  the  inertial  terms  (Equation  Ah)  for  b < h 
as  follows: 

\2 


V*1  . h 
4%  = (» - fl  . 


m 


• 

(1  l\  , P2 

1 f*2 

/l  1 \ P2 

= f 

pi(b  - hj + r 

- 2 1 

V"?)vJ 

(A23) 


b < h (A2l+) 


, b < h (A25) 


7.  In  contrast  to  Figure  A3,  Figure  Ah  illustrates  the  situation 
when  b > h . 


LEGEND 

a = CAVITY  RADIUS 

h = INTERFACE  RADIUS 

b = RADIUS  OF  PLASTIC  FRONT 


Figure  Ah.  Relative  location  of  cav’ by  surface,  layer 
interface,  and  plastic  front  when  b > h 


For  this  case, 


A7 


and 


VV}  * P2 


dr  , b > h 


Evaluating  the  integrals  in  Equations  A26  and  A27  and  substituting  Equa- 
tion A13  into  Equation  A12  result  in  the  inertial  terms  (Equation  A^) 


for  b > h: 


2 (a,b)  = 
■ra  ’ 


aopl  j r V 2 

•*>>’ 


(A28) 


yb>“>  = p2(f-^ir)  >b>h 


, b > h (A29) 


(A30) 


Shear  Terms 


8,  Consider  the  term  I (a,«>)  , which  constitutes  the  effect  of 

S 12 

shear  in  Equation  A4;  the  appropriate  stress-strain  relations  are  ex- 
pressed as 


o.  - a = 
0 r 


r-  E(e_  - e ) , elastic  state 
5 w r 


^ Y + ~ E^Ceg  - c^)  , plastic  -state 


(A31) 


where  E and  E^  are  elastic  and  strain-hardening  moduli,  Y is 
yield  strength,  and  and  eQ  are  radial  and  tangential  strains. 


I 


The  locking  assumption  implies  that- 


li,12 


e •+  2e  = 
r 0 


( 0 , elastic  state 
I , plastic  state 


(A32) 


The  quantities  E , E,  , Y , and  e are  material  constants.  It  is 
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appropriate  to  express  the  tangential  strain  as 


S0  = m (A33) 

o 

where  rQ  and  r are  initial  and  instantaneous  particle  positions, 
respectively.  Combining  Equations  A31-A33  yields 


2E  In 


r > b 


oA  - a = 
0 r 


(A3h) 


Y + 2Et(lnf--ieAj 


r < b 


where  the  elastic  and  locked  plastic  regions  are  separated  by  the  plas- 
tic front  b(t)  . Note  that  the  shear  stress  on  - a does  not  have  to 

Or 

be  continuous  for  Equation  A 6 to  exist;  piecewise  continuity  is  suffi- 

2i 

cient  for  existence.  A discontinuity  in  a0  - ar  exists  at  r = h 
since  this  is  the  material  interface  where  the  material  properties  E , 
E^  , Y , and  are  discontinuous.  Another  discontinuity  in  o0  - 
exists  at  the  plastic  front  r = b(t)  where  the  material  undergoes  a 
sudden  compression.  However,  when  = 0 in  Equation  A3^ , there  is  no 
discontinuity  in  °q  - °r  r = b ^ h . Restricting  attention  to 
situations  where  |e^|  <<  3Y/2E^  reduces  Equation  A3^  to 

2(E  - E ) In  s Y , |ej«f§-  U35) 

o t 


at  r = b(t)  , where  bQ  is  the  initial  position  of  a particle  located 
at  b(t)  . By  rearranging  Equation  A35,  the  quantity  3 is  defined  by 


A9 


-Lil 


»o  " 2<E  - V 


and  it  follows  that 


b sbe 
o 


Conservation  of  mass  in  the  elastic  region  requires  that' 


3 3 3 3 

r ~ b = r - b , r > b 
o o ’ 


Substituting  Equation  A37  into  Equation  A38  and  rearranging  results  in 


r 1 , T b3  -3Bn1 

rc  3 ln£  ~ r3  (1  ‘ 6 ’ 


r > b 


9.  In  the  locked  plastic  region  surrounding  a cavity  of  radius 

11  12 

R(t)  and  initial  radius  R , conservation  of  mass  demands  that  5 


o “ Ro]  = p£(r3  " r3)  » R < r < b 


which  gives  rise  to 


where 


ln  7“  " 3 cl  = ‘ 3 ln  I1  “ ) ’ R < r < b 


R*  = Rll  - ~ 

\ piR 


e = ln  ■£— 
* p* 


Incorporating  Equation  A39  into  the  first  of  Equation  A3H  and  integrating 
with  respect  to  r from  r = rn  to  r = r0  for  constant  E obtains 


3E^~~  , b < T±  < r2  < » (A44) 

with  the  restriction  that  8 « 1 . Likewise,  employing  Equation  A4l  in 
the  second  of  Equation  A34  and  integrating  from  r = r..  to  r = for 
constant  Y and.  E^  results  in 


R 1 rx  1 r2  i 1 (A45) 


(Note:  Equations  A44  and  A45  are  evaluated  in  Reference  11.)  Equa- 

tions A42,  A44,  and  A45  now  make  evaluation  of  Equation  A 6 straight- 
forward. With  a(t)  and  h(t)  representing  cavities  relative  to  the 
first  and  second  layers,  respectively,  R#  , R , and  Rq  are  replaced 
in  the  same  order  in  Equation  A42  by  a#  , a , and  and  then  by 

h-  , h , and  h so  that 
* . ’ o 


where  the  subscripts  1 and  2 denote  first  and  second  layer  properties. 

10.  Examining  first  the  case  where  b(t)  < h(t)  in  Figure  A3 
makes  it  possible  to  rewrite  the  shear  term  in  Equation  A4  as 

I (a,«°)  = I (a,b)  + I (b,h)  + I (h,°»)  , b < h (A4T) 

s s s s 

Applying  Equations  A44  and  A45  to  Equation  A47  then  produces 


All 


CO 

",  ,.3m 

. . 3nT! 

I (a,“) 

5 

* 2h ln  r + 9 \ ^ 

th 

f / m 

(?) 

! 

i 

m=l 


+ I ~$j+3  *1*2  £ ’ " h (M8) 

For  the  case  vhere  b(t)  > h(t)  in  Figure  A4,  Equation  A47  is  replaced 
by 


I (a,-)  = I (a,h)  + I_(h,b)  + I (b,«)  , b > h (A49) 

s s s s 


By  means  of  Equations  A44  and  A45,  Equation  Ak9  is  superseded  by 


+ 3 32E2  , b > h 


(A50) 


Relations  Among  a(t)  , b(t)  , and  h(t) 


11,  First,  examine  Equations  A3?  and  A4o  for  the  case 
b(t)  < h(t)  . Replacing  r and  rQ  by  h and  hp  in  Equation  A38 
and  combining  with  Equation  A37  yields 

' -3B ' 


h3  - li3  = ( 1 - e 1ib3  , b < h 


(A51) 


Likewise,  replacing  r , rQ  , R , and  Pq  by  b , bQ  , a , and  ao  , 


A12 


respectively,  in  Equation  A40  and  combining  vith  Equation  A37  results  in 
the  expression 


a3  " <5.jb3  > b < h 


(A52) 


where 


, n pl  "331 

1 ^7e 

i 


(A53) 


and  a#  is  given  by  the  first  of  Equation  A 46.  Evaluating  Equation  A51 
at  b = h obtains 


h = h e 
o 


b = h 


(A54) 


Now  since  « 1 , it  is  appropriate  to  make  the  approximation 


h a h , b < h 
o’  — o 


(A55) 


By  restricting  attention  to  cavities  which  expand  to  more  than  twice  the 
initial  radius  aQ  , it  is  acceptable  to  replace  the  first  of  Equa- 
tion A46  by 


% ~ a 


(A56) 


so  that  Equation  Af;2  reduces  to 


a3  s , b < h 


(A57) 


a^a  = 6-b^b  , b < h 
1 ’ o 


(A58) 


12,  When  b(t)  > h(t)  , r , rQ  , R , and  Rq  are  replaced 
respectively  in  Equation  A40  by  b , bQ  , h , and  hQ  . Consolidating 


A13 


the  resulting  equation  with  Equation  A37  obtains 


h3  = 6 b3  , b > h 


(A59) 


and 


h2h#  = 6gb^b  , b > hc 


(A60) 


where 


. p2  "362 
62  = ! - ~ e 

2 


(A61) 


Now  Equation  A^O  with  r , rQ  , R , and  R^  replaced  respectively  by 

h , h , a , and  a becomes 
o * o 


pl(ho  - ao)  ■ p«  <»3  - a3>  • 


(A62) 


By  rearranging  and  combining  with  Equations  A46  and  A56,  Equation  A62 
reduces  to 


and 


h3.4  + fl_h3..a3+f^h3 

* Pp  0 Pp  O 


(A63) 


h2h  = h2h#  = Sgb2^  = a2a  , b > h. 


(a64) 


Combined  Effects 

13.  Inserting  Equations  A58  and  A6^  into  Equation  A17  produces 


f(t)  = 


where 


- 2. 
C.^a  a , 

b 

< h 

( 

„ 2. 
0oa  a , 

b 

> h 

t 

(A65) 


Al*l 


C.  = and  = 

1 1 “ “i  2 1 - «2 


Using  Equations  A15,  A58,  and  A65  in  Equations  A23-A25  results  in 

<*,P„  (1  - 6, )2 


v-1  - v,  V» 

• i><ho 


Im(b,<°)  = C^ad  + 2a  ) 


Na  a\  al  1 _2.2 

b " hj  p2  hj  “ 2 Cla 

/ 4 4\  4‘ 

x .pA7"7/  + p27  • 


b < h 


By  means  of  Equation  A56,  Equation  A48  reduces  to 


Is(a}°°) 


b 4 

2Y1  lna+9Et 


Jll  1_  /a' 

, 6 " / 2 Vd, 

1 Z-f  in 


+ I Vif1  - jf)  + I 61E2  J > ”<ho 


where 


vl--: 

2 j.-t 


Incorporating  Equations  A15,  A64,  and  A65  in  Equations  A28-A30  yields 


a2 ^ ^ / 4\ 

A<?  = /a  \ &2  ^ > h 

b «*(!  -a2)  \7y 


m 


r + p« 


/t-t): 


+ p 


/ 4 4' 

/ a a 

*2U1'  ’d\ 


, b > h 


I (b,<»)  = p,. 


m 


n / , 0.2,  a 1 ~2.2  / a 

C0(a&  + 2a  ) ^ - g 'C2a  ( Tj 

b 


, b > h 

* r 


By  utilizing  Equations  A56  and  A59>  Equation  A50  becomes 


Is(a,°°)  '= 


2Y  In  - + £■  E. 

1 a 9 t. 


3m 


m=l 


+ 2Y2  In  - 


+ K 


00 

x 

2 l-J  in 

m=l 


3m 


- 1 


+ 3 32E2  ’ b > ho 


Summary  of  Results 
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14.  For  consistency  with  previous  work,  the  compressive 
stress  p(t)  at  the  cavity  surface  is  expressed  as 


p(t) 


Al6 


(A72) 


(A73) 


(A74) 


(A75) 


normal 

(A76) 


where 


h 


tS^ZSk&Km 


l 


f 


I 


t, 

(P- 


Ps  = Ig(a,») 


and 


,^V*  + B ,/)  = 4<S„  + Im(a,b)  ♦ !„(]>,»> 


For  the  case  where  b.  < h , Ps  > » an<^  ^2  are  ®^ven  ^ 

2 A , ..,3ml 


f"2M» 


m=l 


h' 


a. 


, b < h 


_?B  , “li1  - {i)2  (l\.d. 
^(i-v  Uv  2\ 


1 - 


61  / r> 


2(1  - aj7 


H-Si 


(ATT) 


(A78) 


+ 3W-J  httj’  b<ho  (AT9) 


(A8o) 


, b < h (A8l) 
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where  h^h  and  a/b  - SZ*  from  Equation  A57-  For  the  case  where 
o i 


b > hQ  , ps  , Bn  , and  B£  are  expressed  as 


p = 2Y.,  In  “ + r E, 
*s  1 a 9 t. 
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from  Equation  A63,  and  b is  given  by  Equation  A59*  Note  that  at 
b ='h  - hQ  , p(t)  may  be  discontinuous.  This  arises  from  the  mathe- 
matical treatment  of  the  problem  and  does  not  necessarily  correspond  to 
a physical  discontinuity.  What  happens  is  that  at  b = h , the  speed  at 
which  the  plastic  front  is  propagated  undergoes  a step  change  due  to  the 
discontinuity  in  Y and  E . This  is  an  elastic  effect  and  can  occur 
even  if  the  materials  are  incompressible. 


Special  Cases 

15.  When  both  layers  have  exactly  the  same  properties  (i.e.  the 
homogeneous  case),  then  by  dropping  the  subscript,  Equation  A57  becomes 


a ^ gl/3 


(A86) 


so  that  Equations  A79-A8H  reduce  to 


p3  = - | Y in  « + | Et 
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(A91) 
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15.  When  a = 0 , the  material  is  incompressible , - and  Equa- 
tions A67  and  A68  reduce  to 


B.  = 1 


(A9 


B2  = l 


(A*, 
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Since  it  is  usually  true  that  ‘ Y « E and  E . « E , then 


2 

"s  3 
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10 


in  agreement  with  Goodier.  On  the  other  hand,  when  compressibility 
predominates  over  shear,  3 « o , and 


B±  a 1 - S 
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i 

I and 


B2  “ I “ 6l/3(:L 


+ 6)  + 


1 >/3 


in  agreement  with  Hanagud  and  Ross 
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APPENDIX  B:  COMPRESSIBILITY 


Introduction 


1.  The  cavity  expansion  theory  developed  in  Appendix  A and  the 

penetration  theory  developed  in  Part  II  Loth  utilize  the  concept  of 
11  12 

locking  behavior  * to  estimate  the  effect  of  material  compressibility 
at  any  given  instant.  Moreover,  within  the  context  of  the  cavity  expan- 
sion theory  and  the  present  penetration  theory,  the  incorporation  of 

compressibility  constitutes  a perturbation  of  the  basic  incompressible 

9 10 

theories  set  forth  by  Goodier.  * In  this  appendix  a technique  for 
approximating  the  effective  locking  strain  is  developed  in  which 
the  material  is  first  treated  as  incompressible  in  order  to  obtain  an 
approximate  value  for  the  average  dynamic  pressure  in  the  plastic  zone. 
This  pressure  can  then  be  used  to  estimate  the  average  volumetric  strain 
(i.e.,  the  effective  locking  strain  e&)  in  the  plastic  zone  when  the 
material  is  only  slightly  compressible. 


Average  Dynamic  Pressure  in  the  Plastic  Zone 


2.  Suppose  there  is  a spherical  cavity  centered  at  the  origin 

with  radius  a(t)  , radial  velocity  a , and  radial  acceleration  £ . 

Furthermore,  suppose  that  the  cavity  is  surrounded  by  an  incompressible 

homogeneous  material  whose  behavior  in  shear  is  rate  independent.  The 

situation  is  spherically  symmetric,  and  the  dynamic  pressure  P at  any 

22 

r > a in  the  material  is 


P = 


(aS  + 2a2)  - 


4.2 


a a 


(Bl) 


where  r is  the  radial  coordinate,  p is  density,  and  the  material  is 
stress-free  as  r-*°  . The  average  value  of  P in  a region  a < r < b 


Pr2  dr 


r2  dr. 


Employing  Equation  B1  in  Equation  B2  and  evaluating  the  integrals  re- 
sults in  the  expression 

••• 5 f-fl  . 


Suppose  now  that  b represents  the  radius  of  the  plastic  front  for  an 
incompressible  elastic-plastic  material  as  described  in  Appendix.  A 
(where  p = ) such  that  Equations  A8 6 and  A90  reduce  to 

- dij/3  (Bli 

b i 


= 1 - e 


Equation  Bh  is  now  incorporated  in  Equation  B3  so  the  average  dynamic 
pressure  in  the  plastic  zone  a < r < b for  an  incompressible  material 


Equation  Bb  can  now  be  used  in  conjunction  with  a pressure-density  curve 
to  obtain  a first  approximation  to  the  average  volumetric  strain  (i.e. 
locking  strain)  in  the  plastic  zone  a < r < b as  long  as 

|e^j  « 1 . Thus,  for  materials  which  are  only  slightly  compressible, 
Equation  B6  can  be  used  to  obtain  a value  for  e.  which  is  appropriate 


for  the  cavity  expansion  theory  in  Appendix  A.  For  a two-layer  situa- 
tion, Equation  B 6 is  used  in  turn  for  each  layer.  This  may  lead  to  an 
underestimation  of  |e^|  for  the  first  layer  and  an  overestimation  for 
the  second  layer,  hut  a reasonable  approximation  to  the  average  effect 
of  compressibility  for  both  layers  is  obtained. 


3.  An  upper  limit  on  P&ve  for  projectile  penetration  is  obtained 
by  replacing  a with  the  projectile  velocity  q and  dropping  the  mate- 
rial acceleration  term  a £ (Appendix  C)  in  Equation  B6  such  that 


P-J  «> 

_ O 


/ps1/3  **  + >/3 

2l26I  " 36I  + 6I 


1-6. 


Equation  B7  can  be  used  with  a pressure-density  curve  for  a given  target 
material  to  obtain  an  upper  limit  on  |e^|  at  any  velocity  q as  long 
as  |ej  « 1 . Subject  to  this  limitation. 


a = i - £-  a 


APPENDIX  C:  TARGET  ACCELERATION 


Introduction 


1.  In  Goodier's  original  penetration  theory  for  rigid  spherical 

9 10 

projectiles,  the  projectile  equation  of  motion  is 

(M  + f ’^1)4 - -TOo(*s  + 1 <C1) 

The  development  of  this  equation  is  discussed  in  Part  II,  paragraph  9* 

The  term  (2/3)ta  p0B1  reflects  Goodier’s  assumed  variation  of  target 

O J6  J. 

acceleration  over  the  frontal  surface  of  the  projectile,  and  the  quantity 

*5  TO 

(2/3)ira“p  is  generally  known  as  the  added  mass  for  the  sphere.  When 
0 

R « 1 (Part  II,  Equation  5),  this  term  has  only  a small  effect  on 
£1 

penetration  depth  predictions;  "but  when  the  added  mass  is  comparable  to 

projectile  mass  M (e.g.,  R £ O.l),  then  the  effect  of  this  term 

£1 

becomes  significant. 

The  Effect  of  Added  Mass  on  Penetration  Predictions 

2.  Consider  the  projectile  equation  of  motion  (Cl)  which  without 

O 

the  term  (2/3)tra  p.B.  reduces  to 
o X,  1 

Mil  = -«aB(ps  + § PtB242)  (02) 

where  deletion  of  the  added  mass  is  tantamount  to  disregarding  target  ac- 
celeration. Equation  C2  is  now  used  in  Goodier’s  original  applications 
of  his  rigid  projectile  penetration  theory.  The  results  are  shown  in 
Figures  C1-C4;  the  solid  curves  correspond  to  Equation  Cl,  and  the  dashed 
curves  correspond  to  Equation  C2.  In  these  figures  q^  is  final  pene- 
tration depth,  D is  projectile  diameter,  v is  impact  velocity,  p 

o p 

is  projectile  density,  and  Y is  target  yield  strength.  In  each  case 

discarding  the  added  mass  term  improves  agreement  with  experiment.  From 

3 

these  results  it  is  concluded  that  the  contribution  of  the  terra  ira  p.B, 

o’x  1 


Cl 
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Figure  Cl.  Penetration  of  aluminum  target  (AL  2024  T3) 

by  lead  spheres,  R = 0.l6 
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Figure  C2.  Penetration  of  aluminum  target  (AL  2024  T3) 

by  copper  spheres,  R =0.22 
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C2 


whj  _-h  appears  in  Part  II,  Equation  30  and  elsewhere,  probably  represents 
a high  upper  limit  insofar  as  the  effect  of  target  acceleration  on  final 
depth  of  penetration  is  concerned.  Furthermore,  for  most  applications 
of  the  penetration  theory  as  formulated,  it  appears  appropriate  simply 
to  drop  the  term  from  Equation  30.  However,  this  term  is  retained 
as  an  upper  limit  for  added  mass  throughout  the  presentation  of  theory 
in  Part  II. 


ch 


APPENDIX  D:  NOTATION 

Radial  position,  velocity,  and  acceleration  of  spherical  cavity 
wall,  respectively 

Projectile  radius  at  point  of  contact  with  a layer  interface 

Projectile  radius  at  point  of  contact  with  the  target  surface 

Projectile  radius  at  "base  of  nose 

2 

Projectile  frontal  area  = ira 

o 

Radii  of  spherical  plastic  front  and  concentric  layer  inter- 
face, respectively 

Inertial  coefficients  from  cavity  expansion  theory 

Initial  pound  speed 

Center  of  gravity 

Ogive  caliber  radius  head 

Projectile  diameter  at  base  of  nose  = 2aQ 

Modulus  of  elasticity 

Strain-hardening  modulus 

Factor  relating  projectile  nose  shape  to  target  inertial 
resistance 

Resisting  force  projected  normal  to  projectile  axis  of  symmetry 

Resisting  force  projected  along  projectile  axis  of  symmetry 

Ratio  of  maximum  inertial  resistance  to  shear  resistance  within 
the  context  of  penetration  theory — 

Layer  interface  location 

Initial  layer  interface  location 

Projectile  moment  of  inertia 

Momentum  integral 

General  shear  integral 

Nose  length 

Projectile  mass 

Moment  about  CG 

Compressive  normal  stress  on  projectile  surface 
Average  value  of  p over  the  projectile  nose  surface 
Component  of  p due  to  target  deformation  (shear  resistance) 
Component  of  p due  to  target  inertia  (dynamic  pressure) 


(p£B242/ps^ 
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Average  dynamic  pressure  in  the  plastic  zone 

Projectile  penetration  depth,  velocity,,  and  acceleration, 
respectively 

Final  depth  of  penetration 

Final  penetration,  depth  for  a blunt  cylinder 

Projectile  penetration  depth  when  the  plastic  front  reaches  a 
layer  interface 

O 

i?a^p^/M  , added  mass  parameter 
Nose  performance  ratio 

• P . 

pq  /Y  , order  of  magnitude  index  of  the  ratio  of  inertia  forces 
to  shear  forces,  referred  to  as  the  "solid  Reynolds  number" 

•v 

Location  of  plastic  front 
Time 

Radial  velocity  in  the  elastic  region 
Radial  velocity  in  the  plastic  region 
Normal  component  of  particle' velocity 
Impact  velocity 
Target  particle  velocity 
Tangential  component  of  particle  velocity 
Projectile  weight 
Fixed  space  coordinates 
Yield  strength 

Distance  from  projectile  nose  tip  along  axis  of  symmetry 

Distance  from  projectile  nose  tip  to  CG 

Point  where  the  projectile  axis  of  symmetry  intersects  the 
target  surface 

Target  compressibility  = 1 - p/p^ 

Y/2/(E  - Et) 

Moment  arm 
1 - (p/P£)  exp  (-33) 

Locking  strain  = In  (p/p^) 

Radial  strain 
Tangential  strain 
L/D 

Polar  coordinate  or  polar  angle 


D2 


X 

See  Equation  78 

K 

z/z 

0 

ir 

3.lUl6 

P 

Initial  target  density 

PA 

Locked  target  density,  i.e.-,  average  compressed  density  in  the 
piastic  zone 

4+) 

Limiting  value  of  a<  as  r -*■  b 

r 

in  the  elastic  region 

4-> 

Limiting  value  of  as  r -*■  b 

in  the  locked  plastic  region 

a 

r 

Radial  stress 

V 

Circumferential  stress 

0 

Cone  half-angle 

* 

Circumscribed  cone  half-angle 

a),(o,a)  Angle  of  obliquity,  rotational  velocity,  and  angular  accelera- 
tion, respectively 

oi  Initial  angle  of  obliquity 


Note:  A dot  over  any  quantity  denotes  a derivative  vith  respect  to  time 
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In  accordance  vlth  HR  70-2-3,  paragraph  6c(l)(b), 
dated  13  February  1973,  a facsimile  catalog  c&rd 
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